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Quasi-Optical Polarization Diplexing 
of Microwaves 


By T. S. CHU, M. J. GANS, and W. E. LEGG 
(Manuscript received July 17, 1975) 


The feasibility of a microwave quasi-optical polarization diplexer has 
been demonstrated using photo-etched copper strips with thin mylar backing 
as a practical example. The insertion losses of the principal polarization 
and the cross polarization have been calculated and measured. The con- 
ducting strips must be aligned in a preferred direction, namely, perpen- 
dicular to the plane of incidence, to minimize cross-polarized radiation, 
whereas orientation of the plane of the grid with respect to the beam direc- 
tion rs not restricted. The measured cross-polarized radiation agrees with 
predictions from simple theoretical models of a magnetic current sheet 
for the transmission mode and an electric current sheet for the reflection 
mode. This type of diplexer has been successfully employed in studying 
the polarization properties of the 20-GHz signal from the ATS-6 satellite. 


|. INTRODUCTION 


To achieve frequency reuse by employing orthogonal polarizations 
in a radio communication system, it is essential to avoid cross polari- 
zation in the feed patterns that illuminate the antennas; this relies 
upon the diplexing of two orthogonal polarizations with high isolation. 
Waveguide-type polarization couplers perform diplexing well where a 
carrier with an associated bandwidth of about 10 percent is involved. 
However, it is difficult to provide a low-loss (<0.1 dB) waveguide di- 
plexer to separate effectively and simultaneously the two polarizations 
in each of two widely separated common carrier bands, such as 18 and 
30 GHz. The difficulty stems from the vulnerability of an oversized 
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waveguide to higher-order modes in the higher-frequency band. Con- 
tamination by only 1 percent of the power in higher-order modes may 
cause unacceptable cross-polarized radiation in the feed. 

The necessity for overcoming this problem has led to the suggestion 
of using a closely spaced wire grid as a quasi-optical polarization 
diplexer. The purpose of this paper is to describe feasibility studies of 
this quasi-optical approach. Not only should the insertion loss be 
small and the feed pattern distortion slight in the principal polariza- 
tion, but minimization of the cross-polarized radiation should also 
be achieved. Our measurements have shown that the wires must be 
oriented in a preferred direction to minimize the cross-polarized radia- 
tion. This property is explained theoretically by utilizing a magnetic 
current sheet for transmission through the grid and an electric current 
sheet for reflection from the grid. This preferred direction requires 
that the wires be perpendicular to the plane of incidence determined 
by the beam axis and the grid normal, as shown in Fig. 1b. One notes 
that the classical application of a polarizer, which consists of a wire 
grid parallel to an aperture plane, always satisfies this condition. 

In Section II we calculate the insertion loss and cross-polarized 
radiation using simple theoretical models. Section III describes the 
measurements and the comparison between the calculated and mea- 
sured cross-polarized radiation of wire grids in various configurations. 
Section IV discusses applications and includes concluding remarks. 

To avoid confusion about the definition of cross polarization,’ the 
following two explicit expressions 


é: = 6cos¢ — gsingd (1) 
fp. = sing + dcos¢d (2) 


are defined as the two orthogonal polarization vectors which are the 
cross polarization of each other. The carat ‘‘*’’ indicates unit vector, 
and (6, ¢) are spherical coordinates as shown in Fig. la; #; and fe 
within the immediate vicinity of the Z axis are in the nominal X and 
Y directions, respectively. The usual antenna pattern measurements 
yield directly the patterns for the two polarization components under 
this definition. If a feed pattern with the polarization vector (1) or 
(2) (often called a balanced-feed radiation?) illuminates a paraboloid 
with its axis oriented in the Z direction, the reflected field in the 
aperture of the paraboloid is free of cross polarization. 


ll. THEORETICAL CALCULATIONS 
2.1 Insertion losses 


Implementation of quasi-optical polarization diplexers requires pre- 
diction of the insertion losses for both the principal and the unwanted 
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Fig. 1—(a) Geometry of a quasi-optical polarization diplexer. (b) Preferred con- 
figuration of quasi-optical polarization diplexing. 


cross polarization. If the polarization grid is made of uniformly spaced 
copper strips with a thin mylar backing, the calculation is facilitated 
by the equivalence between the periodic grating and the capacitive 
diaphragm in a parallel-plate waveguide.’ Neglecting the effect of the 
thin mylar layer, the power reflection coefficient for an incident wave 
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polarized perpendicular to the strips is approximately given by? 


B? cos? 6 
4 + B? cos? 6’ (3) 


4b ad 
B= z 2 sec (Sr oT ) 


is the shunt susceptance, b the grating period, d the width of the 
conducting strip, \ the wavelength, and @ the angle of incidence be- 
tween the direction of propagation and the normal to the grating 
plane. Equation (3)* is based upon the low-frequency approximation 
(b<)) for a grating of infinitely thin perfectly conducting strips. 
Using Babinet’s principle, the power transmission coefficient “7,” for 
an incident wave, polarized parallel to the plane determined by the 
strip and the propagation direction, is also approximately given by 
eq. (3) provided the strip width d in the expression for B is replaced 
by the spacing (b — ad). 

If we employ the numerical value, b = 0.56 mm, d = 0.2 mm, 
\ = 1.05 cm, and @ = 45°, as used in the experiment? later, substitu- 
tion into eq. (3) yields the insertion loss of the cross polarization 


—10 logio Ry = 36.8 dB 
=— 10 logio Th 28.8 dB, 


R, = 


where 


which is equivalent to an insertion loss of only 0.001 dB and 0.006 dB 
for the transmitting and reflecting principal polarizations, respectively. 
Typically, the cross polarization in a Cassegrain feed aperture that 
illuminates a polarization diplexer is of the order of —20 dB or less; 
thus, the improvement provided by use of a quasi-optical polarization 
diplexer reduces the residual cross-polarized components in the grid 
aperture to negligible values. 

However, the residual cross polarization discussed above is only part 
of the possible cross-polarized radiation; the following calculations 
show that the co-polarized field in the grid aperture may also give 
rise to off-axis cross-polarized radiation if the strips are not oriented 
in a preferred direction. 


* Equation (3) is known to be accurate when the strips are parallel to the plane 
of incidence for all strip widths. However, to the authors’ knowledge, the rigorous 
demonstration, in the literature, of its accuracy for other strip directions is restricted 
to the cases (b ~ d)/b <1, or d/b «1. 

Equation (8), for the idealized grid, shows that d = 6/2 provides the minimum 
value of the quantity, max [R,, Tu]. Therefore, the grid was designed to have both 
copper strip width and gap spacing equal to 0. 25 mm. However, this specification was 
close to the resolution limit of the fabrication process in use at the time, which re- 
sulted in a grid with the above measured dimensions. 
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2.2 Magnetic current sheet 


For the case of transmission through an arbitrarily oriented wire 
grid, a magnetic-current-sheet model is used to calculate the radiation. 

If a wire grid is placed in front of a radiating aperture that produces 
a wave collimated in the Z direction, the on-axis radiation is linearly 
polarized perpendicular to the conducting direction of the wire grid. 
We shall find the polarization properties of the off-axis radiation. 

Let the plane of the wire grid be oriented in an arbitrary direction, 
as shown in Fig. la, with the following unit normal 


a= #sinBcosa+ ¢7sinBsina + ZcosB. (4) 


In order that an incident electric field, Ey = Hi%, may pass freely 
through the wire grid, the direction of the conducting wires 


ww = siny gy + cos y 2 (5) 


must be the same as that of the equivalent magnetic current density 
2% X Ez (see the appendix) : 


= sin 6 sin a | (6) 
Vceos? 8 + sin? 6 sin? a : 


cos B 
Veos? 6 + sin? B sin? a 


sin y = 





and 
—sin 6 sina 


cos y = : — 
Veos? 6 + sin? B sin? a 








The far-zone electric-field radiation of a magnetic current sheet can 
be written 





je, : : hett ds 
Bi= —srgei® f [RX (x Ey) Jem Aad, (7) 
where & is the free-space phase constant, and R and R are the distance 
to the far-field point and the corresponding unit direction vector. The 
points on the magnetic current sheet are defined by R’. The polariza- 
tion is determined by the bracketed vector product in eq. (7), 





P, = R X (4 X E)) = Ervecos? B + sin? B sin? a [4(—cos ¢ sin y) 
+¢(—sin 6 cos y + cos @sin¢gsiny)]. (8) 
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The dot product of eqs. (1) and (8) gives the principal polarization 
component 


P,-/, = —£, cos B[1 — sin? (1 — cos 6) — sin @sin¢@coty]. (9) 


If we substitute the above product for the bracket in eq. (7), it is seen 
that the on-axis radiation is the same as that without the wire grid, 
while the off-axis radiation is only slightly perturbed. Here the factor 
cos 8 accounts for the larger slanted area of the grid. 

Now the dot product of eqs. (2) and (8) gives the cross-polarization 
component 


P,-f2 = —£1 cos B[sin ¢ cos ¢(1 — cos @) + sin @cos¢ cot y]. (10) 


The cross polarization on axis vanishes, as expected. When the 
direction of the conducting wire is perpendicular to the beam axis, 
l.e., y = 90°, only second-order cross polarization, as represented by 
the first term in eq. (10), is present with maxima in the ¢ = 45° planes. 
This second-order cross polarization is negligibly small for narrow feed 
patterns. However, this residue can become a considerable item for 
broad feed patterns, as will be demonstrated by an experiment de- 
scribed later. 

When the direction of the conducting wire is not perpendicular to 
the beam axis, 1e., y # 90°, the second term in eq. (10) represents 
first-order cross-polarization lobes with maxima in the ¢ = 0 plane. 
This term is kept small if both 6 and (90° — y) are small. Therefore, 
fine adjustment to reduce the residual cross polarization can be 
accomplished by rotation of the grid in its own plane in the case of a 
narrow feed pattern. 


2.3 Electric current sheet 


Next, consider the case of reflection from the wire grid. The reflected 
field is entirely due to the fields radiated by electric currents flowing 
in the wires. If the grid is fine enough, these currents will flow only 
in the direction of the wires. Thus, we may use an electric current 
sheet model to compute the field reflected from the wire grid. To obtain 
perfect reflection from the grid, the direction of the conducting wires 
must be the same as the induced electric current direction ” X Ha, 
where H, is the incident magnetic field for this case. The far-zone 
electric field of an electric current sheet can be written 


ea fiRx (RX @X H)e™-Aad, (1) 
where Zp is the free-space impedance. The polarization is determined 
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by the vector product inside the bracket in eq. (11): 


P, = RX {R X (4 X H:)} = |4 X He |[6(sin 6 cos y 
—cos @sin ¢ sin y) + ¢(—cos¢sin y)]. (12) 


The above polarization is orthogonal to that of eq. (8). The principal 
and cross-polarization components are obtained by the dot products 
of eq. (12) with eqs. (2) and (1), respectively. 


P.- 62 = —Hz2 cos B[1 — sin? 6(1 — cos 0) — sin @sin¢@coty] (13) 
P.- 6, = H2 cos B[sin ¢ cos (1 — cos 8) + sin 6 cos ¢ cot y]. (14) 


Here the relation |# X H.|sin y = A. cos follows the symmetry 
with respect to # between the incident wave and reflected wave of 
which the magnetic field is equal to H2%. 

Since eqs. (13) and (14) are identical to eqs. (9) and (10) except for 
a proportionality constant, the properties of the off-axis cross-polarized 
radiation described in the previous case are also valid for this orthog- 
onal case in reflection. 


il. EXPERIMENT 
3.1 Insertion loss 

The insertion loss of the combinations of a dual-mode horn and a 
wire grid were measured at 28.5 GHz and 19 GHz in both transmission 
and reflection. Figure 2 shows the sketch of the 28.5-GHz experimental 
model. The dual-mode horn has been described elsewhere.> The wire 
grid was made by photo-etching a copper-covered mylar sheet ; copper 
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Fig. 2—Schematic of experimental assembly for 28.5 GHz. 
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strips 0.2 mm wide and 0.018 mm thick are spaced 0.8 mm apart on 
a mylar sheet 0.013 mm thick. The grid, mounted in a circular wooden 
rim, can be rotated in its own plane, which is oriented at 45° with 
respect to the horn aperture. The 19-GHz experiment uses a scaled 
dual-mode horn and 24.7-cm-diameter circular grid with the same 
copper strips and the same mylar sheet described earlier. 

The measured minimum insertion loss on axis, for both transmission 
and reflection, was found to be only about 0.1 dB for the principal 
polarization; the discrepancy with the calculated 0.001 and 0.006 dB 
can be explained by measuring error and slight pattern distortion due 
to diffraction around the grid. 

The maximum insertion loss of the cross-polarized field, on axis, 
at 19 and 28.5 GHz for both transmission and reflection is shown in 
Table I. The symbols || and _L indicate that the grid wires are parallel 
and perpendicular, respectively, to the plane of incidence. 

The measured data are only in qualitative agreement with the 
approximate prediction from eq. (3). However, the effect of the mylar 
sheet (0.013 mm thick with a dielectric constant of 3), imperfect 
polarization of the horn radiation, and diffraction around the grid 
have been neglected in the approximate calculation. It was observed 
that the measured insertion loss of the cross polarization depends 
somewhat upon the spacing between the horn and the grid. 


3.2 Radiation patterns at 28.5 GHz 


The measured cross polarization in the radiation patterns is found 
to be negligible if the conducting strips are aligned in the preferred 
direction normal to the beam. But for the conducting strips in non- 
preferred directions, such as those parallel to the plane of incidence, 
maximum cross-polarized radiation is obtained in the transverse 
planes—AC for transmission and BC for reflection—both perpen- 
dicular to the plane of Fig. 2. To illustrate the predictions of the theo- 
retical models in the preceding section, we present the measured 


Table |— Measured insertion loss of cross-polarized fields 
28.5 GHz 19 GHz 

Grid-Wire 

Position . oar ; 

Transmission Reflection Transmission Reflection 
||* 24 dB 25 32.5 30 
t 28 dB 30 38 34 
Eq. (3) 28.8 36.8 32.3 40.3 


* ; Conducting strips are parallel to the plane of incidence. 
tT Conducting strips are perpendicular to the plane of incidence. 
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transverse plane patterns at 28.5 GHz for four combinations of horn 
polarization and conducting strip directions. 

The transverse plane patterns in Figs. 3 and 4 were measured with 
the radiation transmitted through the grid. In Fig. 3 the horn polari- 
zation is perpendicular to, and the conducting strips parallel to, the 
plane of Fig. 2. The average of the cross-polarization lobe maxima is 
about 20 dB below that of the principal polarization in the same direc- 
tion (@ = 6°), and agrees well with the prediction of eq. (10) [relative 
to eq. (9) with y = 45°] as shown by the dotted curves. In Fig. 4, 


y= 45° 
—10 ¢= 0° 
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Fig. 3—Radiation patterns of a transmitting grid at 28.5 GHz with conducting 
strips parallel te the plane of incidence. 
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Fig. 4—Measured radiation patterns of a transmitting grid at 28.5 GHz with 
conducting strips perpendicular to the plane of incidence. 


the horn polarization is parallel to, and the conducting strips per- 
pendicular to, the plane of Fig. 2. The measured cross polarization 
of less than —40 dB essentially confirms the theoretical prediction of 
negligible cross polarization from eq. (10) (vy = 90°), since the mea- 
suring accuracy of the cross-polarization level is reliable down to about 
—40 dB. 

The transverse plane patterns in Figs. 5 and 6 were measured with 
the radiation reflected from the grid. In Fig. 5 both the horn polari- 
zation and the conducting strips are parallel to the plane of Fig. 2, 
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Fig. 5—Radiation patterns of a reflecting grid at 28.5 GHz with conducting strips 
parallel to the plane of incidence. 


and the measured cross polarization is essentially the same as that 
of the transmitting case in Fig. 3. In Fig. 6, both the horn polarization 
and the conducting strips are perpendicular to the plane of Fig. 2, 
and the measured cross polarization of less than —40 dB is similar 
to that of Fig. 4. Thus, the results show that in employing quasi- 
optical polarization diplexers the off-axis cross-polarized radiation 
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Fig. 6—Measured radiation patterns of a reflecting grid at 28.5 GHz with con- 
ducting strips perpendicular to the plane of incidence. 


can be suppressed only if the conducting wires are perpendicular to 
the beam (i.e., perpendicular to the plane of incidence). 

Owing to limitation of the measuring accuracy, it is difficult to 
measure the second-order cross polarization, the (1 — cos 6) term in 
eqs. (10) and (14), of the wire grid for narrow feed patterns. Therefore, 
we conducted an experiment with a broad feed pattern to check this 
term, which grows rapidly when @ increases. The radiation patterns 
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Fig. 7—Radiation patterns of a small dual-mode circular aperture feed (D/dX = 1.39 
at 16.5 GHz). 


of a small dual-mode horn® (see inset in Fig. 7) were measured with 
and without the wire grid. 

In the absence of the grid, the measured 45° plane patterns of 
principal and cross polarization are as shown in Fig. 7. The measured 
patterns in other planes (not shown) exhibited circular symmetry in 
the co-polarized radiation pattern, and less than —40 dB in cross 
polarization everywhere. When the small dual-mode aperture was 
covered by a wire grid, the measured pattern in co-polarization re- 
mains essentially the same as without the grid; however, the cross 
polarization in the ¢ = 45° plane rises to —26 dB as shown in Fig. 7. 
The calculated cross polarization, which is plotted as a dotted curve, 
shows good agreement with the measured pattern. The cross-polarized 
radiation of a small grid-covered aperture is similar to that of a dipole. 

Thus, the above results demonstrate that the wire grid is a good 
polarizer for large apertures, whereas improper use of the wire grid 
can even enhance the cross-polarized radiation of a small aperture. 
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IV. DISCUSSION 


It has been observed? that there will be no cross polarization in the 
main reflector aperture of an offset near-field Cassegrainian antenna 
provided no cross polarization illuminates the subreflector. This ideal 
condition can be approximately realized by application of a quasi- 
optical polarization diplexer to the feed of an offset Cassegrainian 
antenna with large effective F/D ratio. The quasi-optical polarization 
diplexer can also be used on a symmetrical Cassegrainian antenna as 
demonstrated by its application in an earth-station receiver’:® for the 
20-GHz ATS-6 signal. 

The basic philosophy of the quasi-optical polarization diplexer can 
be simply stated as a cleaning up of the two orthogonal polarizations 
simultaneously just before illuminating the subreflector. This cleaning 
process is especially desirable if the feed is an offset reflector with 
relatively small I'/D ratio. But the conducting wires must be oriented 
in a preferred direction, perpendicular to the plane of incidence, to 
avoid the off-axis cross-polarized radiation. For the broad feed pattern 
of a small dual-mode aperture, the second-order cross-polarized radia- 
tion from a classical polarizer may exceed that of the dual-mode aper- 
ture without the polarizer. The accuracy of the theoretical predictions 
demonstrates the utility of equivalent current sources for such analyses. 

To avoid excessive spill-over loss and pattern distortion, the quasi- 
optical diplexer should be made conservatively large, typically with 
edge illumination less than —20 dB. The main disadvantages of quasi- 
optical feed systems appear to be bulkier volume and heavier weight 
compared with conventional waveguide diplexing feed systems, especi- 
ally when both polarization and frequency diplexing are performed 
by quasi-optical components. But these components have an ad- 
vantage in handling high power without difficulty. Equation (10) in- 
dicates that fine tuning (y 90°) of the residual polarization response 
may be accomplished by rotation of the grid in its own plane. 


APPENDIX 


To calculate the field transmitted through the wire grid, a magnetic 
current equivalent source is chosen, the choice being governed by the 
following considerations. 

Given the tangential electric and/or magnetic field on the bounding 
surface of a source-free region, we may place equivalent sources on the 
bounding surface to correctly reproduce the original field in the source- 
free region :° 


(¢) Magnetic current: K = E; X #, backed by a perfect electric 
conductor on the bounding surface. 
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(it) Electric current: J = % X Hz, backed by a perfect magnetic 
conductor on the bounding surface. 
(zit) Combination: 


K = E; X #, magnetic current ae 
operating in 


J = % X H,, electric current free space. 


All three equivalent sources give identical results if they are based 
on the true fields, E, and H;. However, in the case of transmission 
through a grid, we do not know the true magnetic field H; on the 
source-free side of the grid. If the polarizer is fine enough, one can be 
sure, though, that the tangential electric field is perpendicular to the 
wires. Thus, the field transmitted through the grid can be predicted 
most accurately by the magnetic current equivalent source backed by 
an electric conducting plane on the grid. Since a tangential magnetic 
current imaged in an electric conducting plane is equal to itself, we 
may include the effect of the electric conducting plane by using twice 
the magnetic current, 2K, operating in free space. 

The far field radiated by the magnetic current density, 2K, in free 
space is 
_ _ ik 


E; ee f LOK) x Rye Aaa, (15) 
gr. 


gaps 


If the grid spacing is very small compared with wavelength, then, 
as the magnetic current 2K = 2E, X # varies from zero on the grid 
wires to maximum in the space between, the other terms in the inte- 
grand of eq. (15) are essentially constant. Thus, we may replace the 
fluctuating K with its average value, Kay,, 


B, = — eat f [0K ug) X Rem Add. (16) 
4rR 
Although we know the direction of K, we do not know its magnitude 
unless the reflection coefficient of the grid is known. In the usual case, 
the grid is designed to introduce negligible insertion loss for the de- 
sired polarization, whence E; on axis should equal that present when 
no grid is used. In this case, the magnitude of 2K,,, would have to be 
such that 
2Kavg = 2E; X # (negligible insertion loss) (17) 


in order that eq. (16) will result in the correct on-axis value for E1. 


By substituting eq. (17) into eq. (16), we arrive at eq. (7), the desired 
equation for computing the field transmitted through the wire grid. 
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Material Structure of Germanium-Doped 
Optical Fibers and Preforms 


By H. M. PRESBY, R. D. STANDLEY, J. B. MacCHESNEY, 
and P. B. O’CONNOR 


(Manuscript received July 8, 1975) 


The structural characteristics of preforms and optical fibers fabricated 
by modified chemical vapor deposition were studied by optical, interfer- 
ence, and scanning electron microscopy. It was observed that the structural 
features resulting from the deposition process are preserved through subse- 
quent processing and appear in the fiber with the exception of a region at 
the center of the fiber. Here, selective evaporation of dopant material from 
the onner surface of the deposit results in a refractive index depression on 
the axis of the optical waveguide. 


I. INTRODUCTION 


The chemical vapor deposition process, in which oxides are de- 
posited and simultaneously fused on the inner surface of a fused silica 
tube, has become a valuable technique for fabricating low-loss!? and 
graded-index optical fibers.’ In modifications and refinements of this 
technique, higher depositional rates and very low-loss single-mode 
fibers* have also been achieved. 

An important question that arises in utilizing this process concerns 
the correlation of the deposited material structure in the preform to 
that in the resulting optical fiber. Can one be confident, for example, 
that the same distribution of refractive index that is introduced into 
the preform by changing the material composition of the deposited 
layers exists in the fiber pulled from this preform? This determination 
is necessary if one is to reliably fabricate those graded-index profiles 
required to achieve a minimum of pulse dispersion.® This is due to the 
fact that the shaping of the index profile is quite critical because the 
reduction-in-pulse-dispersion-vs-profile curves exhibit a singularity- 
like behavior in the region of the optimum index distribution.® 

Evidence for the preservation of the deposited profile has recently 
been reported, based on the observation of a linear increase in refrac- 
tive index in a fiber which was pulled from a preform in which the 
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dopant concentration was increased in the same manner.’ In this 
paper, we present results of optical, interference, and scanning electron 
microscope studies’ of a graded, near-parabolic, index fiber and pre- 
form as further aid in understanding the transition of material from 
the preform to the fiber state. A main conclusion of this study is that 
the structural features resulting from the preform deposition process 
are preserved and, after suitable scale transformation, appear in the 
fiber. Due consideration should therefore be given to depositional 
characteristics that may ultimately affect transmission behavior. 


Il. OBSERVATIONS 


The preform originates from an approximately 4-m-long 12 X 14- 
mm fused-quartz tube which is collapsed into a rod after the deposition 
process. In the structure studied here, the deposition started with an 
initial layer of borosilicate to prevent impurity diffusion into the core. 
This layer was deposited with 41 traversals of a oxyhydrogen burner 
which provides the heat to react the BCl;, SiCl., and Oe starting 
materials. The core deposition process consists of systematically in- 
creasing the flow of GeCl, while holding the flows of SiCl, and BCI; 
constant, thus producing an increasing GeO, content and associated 
increased refractive index with increasing deposit thickness. The 
GeCl, flow was increased 11 times, in such a manner as to produce a 
near-parabolic index variation from the cladding interface to the 
center of the core. The number of torch traversals during each of the 
11 steps was controlled to make the thickness of each step approxi- 
mately equal. After collapse, a length of preform was pulled into a 
fiber with an overall diameter of ~100 um by the use of an electric 
furnace. 

A slice transverse to the axis of the remaining length of preform was 
made and then polished to a thickness of approximately 10 mills for 
interference and optical microscopic observations. Transverse samples 
of the fiber were also prepared for interference-microscope and scan- 
ning-electron-microscope studies.’ In the latter case, after a short 
length of fiber is scored and broken to ensure a flat end, it is etched in 
a 25-percent solution of hydrofluoric acid for several minutes and then 
flash-coated. The last step is performed to prevent charge build-up on 
the sample during scanning-electron microscope observations. 

An overall view of the preform sample observed with conventional 
optical microscopy is shown in Fig. la. The sample is ~7.2 mm in 
diameter with a core diameter of 4 mm. The irregular shape is due to 
pieces of the cladding which broke off during the cutting and polishing 
procedure. It should be noted that the preform is under considerable 
stress because of the difference in the expansion coefficients of the 
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Fig. 1—Transmitted light photomicrographs of transverse section 
section at higher magnification. 


of preform. 
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BORON CLADDING 
LAYER 


(a) Entire section. (b) Central core region. (c) Entire 


core and the cladding. If this stress is not to be relieved by shattering 
upon cutting, extreme care must be used. In this case, the preform was 
cut with a diamond wire saw with controlled lubrication. Even so, the 
preform did crack, as seen by the dark curved line in the left-hand 
section of the sample. 

Figures 1b and 1c obtained by optical microscopy show expanded 
views of the central region and of a portion of the entire cross section. 
The cladding, the borosilicate layer, and the first five germania- 
borosilicate steps are labeled. The regions between the steps are quite 
distinct, as are the individual layers within each step. Each of these 
layers, as noted, corresponds to a traversal of the oxyhydrogen burner 
along the tube. During step 3, for example, the GeCl, flow was main- 
tained constant for nine traversals of the burner, producing the nine 
layers observed in Fig. Ic. 

An expanded view of one-half the preform sample as observed by 
interference microscopy is shown in Fig. 2. The refractive-index differ- 
ence between the cladding and a point in the core of the preform is 
given by the fringe displacement at that point times the wavelength 
of observation and divided by the thickness of the sample. One ob- 
serves the straight parallel fringes in the fused silica cladding on the 
right indicating the uniform composition of this region, as expected. 
The drop in the level of the fringes indicates the termination of the 
cladding and the start of the borosilicate step that has a lower index 
of refraction than pure fused silica. Again, in this region, which extends 
for about 375 um, the composition is relatively uniform and no evi- 
dence is seen of the 41 layers which comprise this step. This tends to 
indicate that some boron diffusion occurs, smoothing out the indi- 
vidual layers. It does not appear, however, that boron diffuses into 
the cladding, as evidenced by the relatively sharp transition occurring 
~10 um between the cladding and this step. 

At the termination of the pure borosilicate layer, germania deposi- 
tion commences. The first three of these steps are labeled. Note, in 
particular, that step 3 exhibits nine sinusoid-like variations which, as 
discussed previously, correspond to nine torch traversals used in 
depositing this step. The reason for this index variation within each 
layer may be due either to a difference in composition of the particles 
reacting homogeneously (i.e., in the gas phase) and heterogeneously 
(i.e., on the surface of the tube),’ or to temperature variation effects, 
depositing different concentrations as the torch passes. These com- 
positional variations become quite sizable as the number of layers in 
each step decreases towards the center of the core. Note that in all 
steps the integrity of these layers is maintained and that relatively 
sharp boundaries exist between the steps. The transitions are most 
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clearly seen between the first few steps in which the change in ger- 
manium content is largest. These observations indicate that little, if 
any, germanium diffusion occurs between adjacent steps or within a 
given step itself. 

The situation is quite different at the center of the core. Despite the 
fact that the germanium concentration was varied in a smoothly in- 
creasing manner in the last several steps, a large disturbance in the 
center is observed, with a corresponding dip in the refractive index. 
We believe this is due to the evaporation of germania at the inner 
surface of the deposit during the elevated temperatures experienced 
in the collapse process. The fact that the last two or three layers 
appear to broaden towards the center also indicates the existence of 
some germanium flow extending beyond the immediate central seg- 
ment. It may be possible to compensate for this effect by a germanium 
overdoping in this region. 

Quantitative index measurements indicate a maximum index differ- 
ence between core and cladding of approximately An,, = 0.016 and an 
index difference between the cladding and the borosilicate layer of 
An = 0.004. Both values are in very good agreement with the corre- 
sponding measurements made from the microinterferrogram of the 
fiber, which is shown in Fig. 3a, and indicate a preservation, for the 
most part of the material composition through the pulling process. 

Note in Fig. 3a the uniform cladding containing straight parallel 
fringes, the subsequent drop owing to the borosilicate layer, and the 
gradual grading of the index profile to a maximum near the axis. 
The grading was near-parabolic and has been related to a reduction in 
pulse dispersion for this fiber.2 These regions are again shown in the 
transmitted-light photomicrograph of Fig. 3b. The core which appears 
as the bright central area is surrounded by the relatively dark borosili- 
cate layer. Beyond that is the grayish cladding with an overall diameter 
of 96 um. 

The resolution of the interference microscope is not sufficient to 
resolve any layer or step structure in the fiber. To obtain greater reso- 
lution, we made use of scanning electron microscopy. Preparation of 
the fiber samples was described previously, and results of observations 
are shown in Figs. 4 and 5. 

Figure 4 presents three micrographs taken at increasing magnifica- 
tions centered on the axis, and Fig. 5 is a composite photograph of a 
section of one-half the fiber at somewhat greater magnification. The 
main points to be noted are the preservation of the step and layer 
structure in the fiber and the appearance of the elevated region near 
the axis. This feature is a region that did not etch as rapidly as the 
surrounding area, because of a lack of contained dopant and agrees 
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Fig. 4—Scanning-electron-beam-microscope photographs of fiber at increasing 
magnification. 


with the preform observation of germanium departure during collapse. 
The appearance of the distinct step and layer structure displays the 
further lack of germanium diffusion even during the pulling process 
and indicates the maintenance of compositional and structurai in- 
tegrity from the deposition through the pulling process, with the ex- 
ception of the central region. 

We further investigated this latter region in another fiber prepared 
in the same manner by chemical vapor deposition. In this fiber, the 
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Fig. 5—Composite scanning-electron-beam-microscope photograph of segment of one-half fiber. 
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GeCl, flow was increased at each of ten steps to produce a linear varia- 
tion in concentration. The number of layers in each step varied some- 
what from the fiber considered above.’ Scanning-electron-microscope 
photographs of etched fiber samples are shown at increasing mag- 
nification in Fig. 6. 

The depression of the index on the axis because of the loss of dopant 
is quite pronounced in this fiber and appears as the micron-or-so in 





(c) 41am 


Fig. 6—Scanning-electron-beam-microscope photograph of fiber having linear 
refractive-index profile. 
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diameter, raised, tapered pip in the center. Note also the preservation 
of the distinct step and layer structure. The refractive index depression 
is also observed in this case in the microinterferrogram of the fiber 
shown in Fig. 7. It is seen as a dip in the fringe normally passing 
through the center of the fiber. Slight modulations of the refractive 
index can also be observed at each of the ten steps in this case owing 
to the relatively large change in germanium concentration between 
the steps. 

It is important to note that, despite the structural features that 
exist on such a small scale, the losses of these fibers were less than 5 
dB/km in the region of 1.0 um. This is presumably due to the fact 
that these features are very uniform in the direction parallel to the 
axis of the fiber and hence do not contribute in a large way to scatter- 
ing losses. Small-scale variations of such features if existing, however, 
could form a lower limit on losses achievable with fibers fabricated 
by this technique. 





[-——— 20 um 


Fig. 7—Microinterferrogram of linear refractive-index profile fiber showing index 
depression in center. 
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In summary, structural features resulting from the preparation of 
both the preforms and graded-index optical fiber by the chemical 
vapor deposition process have been observed by optical, interference, 
and scanning electron microscopy. These features can be directly 
related to steps in the fiber fabrication. It was observed that structures 
present in the preform were preserved through the drawing and were 
present in the fiber. For instance, optical and interference observations 
indicated that germania concentrations varied within each deposited 
layer, and this variation was not substantially altered by subsequent 
processing steps except in the centermost layers. Here, flow of the 
deposit during collapse and vaporization of germania, probably as 
GeO, led to a depleted region of lower refractive index at the fiber 
center. 
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An Efficient Linear-Prediction Vocoder 


By M. R. SAMBUR 
(Manuscript received April 29, 1975) 


A primary interest in any method for producing synthetic speech ts to 
minimize the number of bits per second required to generate acceptable 
quality speech. An efficient method for transmitting the linear-prediction 
parameters has been found by using the techniques of differential PcM. 
Using this technique, speech transmission is achieved employing fewer 
than 1500 bits/s. Further reductions in the linear-prediction storage re- 
quirements can be realized at a cost of higher system complexity by trans- 
mission of the most significant eigenvectors of the parameters. This tech- 
nique in combination with differential PcM can lower the storage to 
1000 bits/s. 


I. INTRODUCTION 


The method of linear prediction has proved quite popular and suc- 
cessful for use in speech compression systems.!~ In this method, speech 
is modeled as the output of an all-pole filter H(z) that is excited by a 
sequence of pulses separated by the pitch period for voiced sounds, or 
pseudo-random noise for unvoiced sounds. These assumptions imply 
that within a frame of speech the output speech sequence is given by 


s(n) = a ays(n — k) + Un, 


where p is the number of modeled poles, wu, is the appropriate input 
excitation, and the a,’s are the coefficients characterizing the filter 
(linear prediction coefficients). Figure 1 illustrates the frequency- 
domain, as well as the equivalent time-domain, model of linear-pre- 
diction speech production. To account for the nonstationary character 
of the speech waveform, the parameters a, of the modeled filter are 
periodically updated during successive speech frames.* Generation of 
speech in this method requires a knowledge of the pitch, the filter 


“A frame is a segment of speech thought adequate to assume stationarity of the 
speech process. Typical frame lengths employed range from 10 to 30 ms. 
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Fig. 1—Discrete model of speech production as employed in linear prediction. 


parameters, and the gain of the filter (amplitude of input excitation) 
in each speech frame. 

A primary interest in any method for producing synthetic speech 
is to minimize the number of bits per second needed to generate ac- 
ceptable quality speech. The smaller the information storage require- 
ments (bits per second), the more attractive the system becomes for 
the important applications of voice answer-back and speech trans- 
mission.’ To achieve the minimum storage requirement for a given 
system, an efficient means of quantizing the generating parameters 
must be determined. Using conventional pulse code modulation (pcm) 
techniques in which the amplitude of each parameter is uniformly 
quantized into 27 levels, it has been found necessary to allot at least 
five bits (B = 5) of information for both pitch and gain and at least 
11 bits for each a,.! The corresponding storage requirements for this 
method of quantization of the linear-prediction (LPc) parameters is 
unacceptable for many applications, and an improved scheme for 
quantizing the parameters is needed. 

For the usual 12-pole linear-prediction representation, the dominant 
portion of storage is allotted to the filter coefficients (132 bits per 
frame in the pcm method of information transmission). The extremely 
fine quantization of the a,’s is necessary because small perturbations 
in the coefficients can sometimes cause radical changes in the important 
pole frequencies of the modeled filter H(z) and may even cause the 
filter to become unstable (poles outside the unit circle). To overcome 
the limitations of quantizing the predictor (filter) coefficients, it has 
been found quite profitable to transmit different but informationally 
equivalent sets of parameters.*:® The most suitable parameters have 
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been experimentally determined to be the log-area ratio coefficients 
g:i.4 These coefficients are nonlinearly related to the filter coefficients by 


L+hk; 
g: = logge, (1) 





where the k,’s are termed the parcor coefficients.!:2-4. If we denote a{” 
as the ith linear prediction coefficient for a jth-pole linear-prediction 
model, then 


k; = af’. (2) 
The parcor coefficients have the very important property that if 
[ke | md 1=1, ot yD, (3) 


then it is guaranteed that the linear-prediction filter is stable.* Thus, 
small perturbations in the parcor coefficients or log-area coefficients 
will not affect the stability of the modeled filter. 

Since the log-area coefficients are nonlinearly related to the filter 
coefficients, transmission of the log-area parameters is equivalent to a 
nonuniform quantization of the linear-prediction coefficients. By 
transmitting the log-area parameters, the number of bits allotted to 
the filter parameters can be effectively reduced by nearly $.3:4-6 In this 
paper, we offer two additional methods of quantization of the necessary 
synthesis parameters (pitch, gain, and filter coefficients) that can even 
further reduce the storage requirements of a linear-prediction vocoder. 
One proposed method of quantization uses the technique of differ- 
ential pcm (ppcm) to transmit the linear-prediction parameters. This 
scheme exploits the fact that the Lpc parameters are themselves pre- 
dictable from previous samples. Using this method, speech transmis- 
sion that is practically equivalent to the unquantized synthesis can 
be achieved using fewer than 2000 bits/s. 

The second method of transmission exploits the redundancy between 
the linear-prediction parameters. The Lpc parameters can be predicted 
not only from the given parameter’s past values, but also in some sense 
from values of the other parameters. The suggested scheme involves 
the transmission (using pPcm techniques) of the most significant eigen- 
vectors of the log-area parameters. For the typical speech utterance, 
the space of the 12 log-area coefficients can be effectively represented 
by only the first five or six eigenvectors. The transmission requirement 
for this method is fewer than 1200 bits/s. 

The organization of this paper is as follows. In Section II, we briefly 
discuss the concept of ppcm coding. In Section III, we show that ppcm 
coding offers a significant improvement over pcm coding for trans- 
mission of the linear-prediction parameters. In Section IV, the results 
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are presented of a synthetic speech experiment using the proposed 
ppcm scheme. In Section V, we discuss several methods of Dpcm coding 
that are more suitable for real-time implementation. Included in this 
section is a discussion of adaptive quantization (appcm) and adaptive 
ppc prediction. In Section VI, we discuss the method of orthogonal 
linear prediction. The results of synthetic speech experiments are 
included in this section. Finally, we conclude with a summary and 
discussion of the results presented in the paper. 


i]. DIFFERENTIAL PULSE CODE MODULATION 


The idea of differential pcm is similar in philosophy to the concept 
employed in linear-prediction speech analysis. In ppcm, we assume that 
the transmitted parameter in a given frame of interest can be estimated 
by a linear combination of the parameter in previous frames.’ If we 
let x, denote the value of the transmission parameter x in the rth 
frame (where x can represent pitch, gain, log-area coefficients, or what- 
ever), then this assumption implies 


tr Vt, = > Bitr—4, (4) 
where 7 is the order of the prcm prediction analysis. The ppcm tech- 
nique calls for the transmission of the difference between the predicted 
value £, and the true value ,. 

Figure 2 illustrates the structure of the ppcm coding system. In the 
implementation of a ppcm scheme, a feedback path around the quan- 
tizer is used to ensure that the error in the reconstructed (quantized) 
signal %, is precisely the quantization error for the difference signal 
€, = X, — #,, where Z is the predicted value based upon the quantized 
signal ,. The predictor coefficients b; are chosen to minimize the power 
of the difference signal e,. The mathematical analysis required for the 
solution of the optimum set of 6,’s is exactly the same as the analysis 
for the calculation of the linear-prediction coefficients, a;,2 = 1, «++, p. 
The determination of the b;’s is made by solving the familiar correla- 
tion equations: 


n N N 
pS b; = Cp_iLr_k FT oe LrLr—ky 1 < k 
r=n r=n 


t=1 


IIA 


n, (5) 


where JN is the number of frames in the utterance. 

The advantage of ppcm coding is obvious when one realizes that, 
if x, can be accurately estimated from previous samples, the informa- 
tion necessary for transmission (as expressed by the difference signal 
x, — £) is necessarily less than the information required for coding 
the signal without exploiting its predictability. The advantage of 
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Fig. 2—Differential pom (az = quantizer; pr = predictor; = 2, = Yi4-, bjty_;). 


DPCM coding can be precisely specified by noting that, for a given fine- 
ness of quantization, the quantization error is proportional io the 
variance of the signal present at the quantizer.’ Thus, the improve- 
ment in performance (as measured by the frequently used standard 
of signal-to-quantization-error ratio) using ppc strategy over straight 
PCM coding is given by the ratio of the variance (power) of 2, to that 
of the difference signal 
_ __ (a7) 

o= ,— By o 
Using the optimum predictors b;, the resulting gain over PcM is ap- 
proximately* given by 





_ 2.9% bC;\7 _ (a) 
Gow = (1 BE) ag ae 7) 
where 
N 
Cc; = ‘y UrLr—i- (8) 


r=n 


For equal standards of synthetic speech quality, the gain obtained 
by using a ppc strategy over that of pcm coding can be traded off to 
reduce the rate of information transmission. Of course, for G < 1, it 
is disadvantageous to use ppcm coding. However, for the transmission 
of parameters that are reasonably smooth in their variation from one 
transmission frame to the next, it is guaranteed that prcm coding is 
superior to Pom coding. In the next section, we demonstrate the effi- 
ciency of ppcm techniques for the coding of the linear-prediction speech 
parameters. 


fil. DPCM IMPROVEMENT IN CODING LPC PARAMETER 

To illustrate the efficiency of ppcm techniques in the coding of the 
synthesis parameters, Fig. 3 shows the improvement factor Gop in 
decibels as a function of the number of ppcm predictors. The figure 


“The gain is approximate because the effects of the quantizer in Fig. 2 are ignored. 
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Fig. 3—Gopt for the sentence, ‘‘“May we all learn a yellow lion roar.”’ 


shows Gop: for the first two log-area coefficients (g: and ge),* pitch 
period and power’ for the all-voiced utterance, ‘‘May we all learn a 
yellow lion roar.’”? The improvement factor was calculated by con- 
sidering each particular parameter across the entire sentence and then 
calculating the optimum predictors using eq. (5) and Gop, using eq. 
(7). The results depicted in Fig. 3 are for a male speaker, but the results 
are typical of those obtained for other male and female speakers. For 
the complete ensemble of parameters necessary to produce synthetic 
speech (12 log-area coefficients, pitch, and power),* the set of improve- 
ment factors were all significantly greater than 1. 

Figure 4 shows a typical plot of the improvement factor calculated 
for a sentence containing unvoiced sounds, “‘Few thieves are never 


*The parameters were calculated at the rate of 50 samples per second. The filter 
parameter was calculated by the covariance method (Ref. 1), and pitch was measured 
by a method developed by B.S. Atal (Ref. 8). 

t Power is defined as the energy in the speech frame. For the synthetic system 
employed, it is more convenient to transmit power instead of the amplitude of the 
input excitation. 

t Log-area coefficients were transmitted because of their optimum quantization 
properties. 
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Fig. 4—Gopt for the sentence, “Few thieves are never sent to the jug.” 


sent to the jug.” In this sentence, the ppcm improvement over PcM 
coding is not as dramatic as for the all-voiced sentence. The reason 
for the decreased values of Gop; is that the synthesis parameters tend 
to change sharply during the unvoiced-voiced transition. Thus, during 
the transition region there is an abrupt reduction in the correlation 
between successive samples, and very little information can be gained 
about the signal from past values. Another reason for the reduced 
values of Gop: 1s that the variation of the Lpc parameters during un- 
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Fig. 5—Gopt for the sentence, “Few thieves are never sent to the jug.”’ A separate 
DPCM analysis is used in each unvoiced and voiced segment. 
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voiced sounds is inherently more random than during voiced sounds 
and is thus less predictable. Fortunately, during unvoiced regions the 
quality of the synthesized speech is more tolerant to quantization 
noise than during voiced regions. Thus, the diminished values of the 
Gopt’S 1S not as significant as might at first appear. 

One method of increasing the improvement factor for utterances 
containing unvoiced sounds is to update the prcm predictors when- 
ever the spectral properties of the speech signal change from unvoiced 
to voicec sounds. Figure 5 shows Gop: for the same sentences as were 
used to obtain the results of Fig. 4, but in this figure the optimum 
ppcm predictors were separately calculated for each different section 
of unvoiced and voiced speech. The improvement factor for this form 
of ppcM coding is about 5 dB better than a single calculation of the 
predictors. In a later section of the paper, we discuss another method 
for updating or adapting the prc predictors to the changing spectral 
properties of the speech signal. 


lV. SPEECH SYNTHESIS 
4.1 Synthesizer 


The improvement factors for the Lpc parameters demonstrate that 
DPcM coding is superior to Pcm coding. To confirm the results of the 
Gopt experiments, a synthetic speech system was constructed in the 
manner illustrated in Fig. 2. To take advantage of the fact that the 
improvement factor saturates near n = 1 (Figs. 3 and 5), only a simple 
first-order ppcm system was used. The optimum predictor was recom- 
puted for each separate unvoiced and voiced region and the LPc 
parameters were calculated at a rate of 50 samples per second. The 
speech was synthesized using the formulation discussed by Atal and 
Hanauer.! After quantization, the parameters were geometrically 
interpolated (linear interpolation on a logarithmic scale) to allow 
pitch synchronous resetting of the synthesizer. 

The quantizer used in the ppcm coding of the synthesis parameter 
was a nonuniform quantizer that was designed to exploit the properties 
of each parameter’s error signal. An experimental investigation has 
indicated that the difference signal for pitch, power, and g; are most 
suitably modeled by a zero mean gamma density, 


P.,(€r) = spe (—kle,|), (9) 


where 
0.75 


CSS * 


k 
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The higher-order log-area coefficients are more Laplacian in character: 


pate) = dom (~ Hel), 


o = v26. 


where 


A signal with a gamma distribution is highly concentrated near its 
mean, but can also readily achieve values more than three standard 
deviations from its mean. A Laplacian signal is less concentrated than 
a gamma signal near its mean value. Figure 6 illustrates the statistical 
characteristics of a zero mean, unit standard, deviation signal with a 
gamma density, a Laplacian density, and a gaussian density. Figure 
7 shows a comparison between the calculated distributions for the 
difference signal of several typical synthesis parameters and their 
approximated distributions. 

For a gamma-behaved signal, the properties of the optimum quan- 
tizer are summarized in Table I.° The x; values in the table define the 
ends of quantizer input ranges, and the y; values are the corresponding 
outputs. Thus, for a two-bit quantizer, an input between 0 and 1.205 is 
quantized as 0.302. Similarly, an input between 0.229 and 0.588 for a 
four-bit scheme is quantized as 0.386. The properties of the optimum 
quantizer for Laplacian signals are summarized in Table IT.9 Included 
in these tables is the expected mean square between the difference 
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Fig. 6—Comparison of a gaussian, gamma, and Laplacian density with zero mean 
and unit standard deviation. 
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signal and the quantized difference. Thus, for a four-bit quantization 


of a gamma signal, the mean square error is 0.0196. 


Tables I and II are constructed for signals with unit standard devia- 
tion. To obtain the levels y; and boundaries «; for signals with standard 

















Table |— Optimum quantizers for signals with gamma density 
(u = 0, a”? = 1) 

B 1 S 3 5 

t Xi Yi Xi Yi Xi Yi Li Yi Li Yi 

1 «© 0.577 | 1.205 0.302 | 0.504 0.149 | 0.229 0.072 | 0.101 0.033 
2 co 2,108 | 1.401 0.859 | 0.588 0.386 | 0.252 0.169 
3 2.872 1.944 | 1.045 0.791 | 0.429 0.334 
4 «© 3.779 | 1.623 1.300 | 0.630 0.523 
5 2.372 1.945 | 0.857 0.737 
6 3.407 2.798 | 1.111 0.976 
7 5.050 4.015 | 1.397 1.245 
8 «© 6.085 | 1.720 1.548 
9 2.089 1.892 
10 2.517 2.287 
11 3.022 2.747 
12 3.633 3.296 
13 4.404 3.970 
14 5.444 4,838 
15 7.046 6.050 
16 2© 8.043 
MSE 0.6680 0.2326 0.0712 0.0196 0.0052 
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Table I] — Optimum quantizers for signals with Laplace density 














(u = 0, o? = 1) 

B 1 2 3 4 5) 

t xi Yi xi Yi xi Yi Xi Yi Li Yi 

1 oo 0.707 | 1.102 0.3895 | 0.504 0.222 | 0.266 0.126 ; 0.147 0.072 
2 eo 1.810 | 1.181 0.785 | 0.566 0.407 | 0.802 0.222 
3 2.285 1.576 | 0.910 0.726 | 0.467 0.382 
4 ) 2.994 | 1.3817 1.095 | 0.642 0.551 
5 1.821 1.540 | 0.829 0.732 
6 2.499 2.103 | 1.081 0.926 
7 3.605 2.895 | 1.250 1.136 
8 20 4.316 | 1.490 1.3865 
9 1.756 1.616 
10 2.055 1.896 
11 2.398 2.214 
12 2.804 2.583 
13 3.305 3.025 
14 3.978 3.586 
15 5.069 4.371 
16 oo 5.768 
MSE 0.5 0.1765 0.0548 0.0154 0.00414 


deviation different from unity, simply multiply the given values by 
the actual standard deviation.* The standard deviation for each param- 
eter can be approximated as the rms power of the unquantized error 
signal. The rms value of the unquantized error signal is obtained 
directly from the calculation of the optimum prc predictors and is 
given by 


4.2 Experimental results 


Four sentences were synthesized in the experimentation: 


A. Few thieves are never sent to the jug. 

B. May we all learn a yellow lion roar. 

C. It’s time we rounded up that herd of Asian cattle. 
D. Should we chase those young outlaw cowboys? 


High-quality recordings of these sentences were made by two male and 
two female speakers, and these utterances were used to obtain the 
analysis data for the ppcm coding method. 


*To obtain the mean square error, multiply the values by the signal variance. 

} Since the properties of the unquantized error signal are explicitly known, it is 
sometimes advantageous to use a more complex nonuniform quantizer to truly 
optimize the transmission system. 
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Various schemes were tested for assigning bit rates for each indi- 
vidual error signal. From informal listening experiments, it was de- 
termined that synthetic speech that was negligibly different from the 
unquantized synthesis could be generated according to the following 
bit assignment : 


Pitch: 3 bits/frame 
Power: 3 bits/frame 
Unvoiced-voiced: 1 bit/frame 
gi: 4 bits/frame 

go: 4 bits/frame 

gs: 4 bits/frame 

gs: 4 bits/frame 

gs: 3 bits/frame 

gs: 3 bits/frame 

gz: 2 bits/frame 

gs: 2 bits/frame 

go: 2 bits/frame 

gio: 1 bit/frame 
gui: 1 bit/frame 
giz: 1 bit/frame 


The total number of bits dedicated to the complete set of LPc param- 
eter is only 38 bits/frame or 1900 bits/s. On the average, an additional 
100 bits/s are required to transmit the necessary ppcm information 
(bPcm predictors, standard deviations, and initial values of the Lrc 
parameters). As can be observed from Figs. 8, 9, and 10, the spectro- 
gram of the ppcm synthetic speech closely resembles that of the un- 
quantized synthetic speech but requires only a fraction of the storage. 

As the bit rate for the ppcm linear prediction vocoder is lowered 
below the value of 2000 bits/s, the quality of the synthesis slowly 
begiris to deviate from that of the unquantized synthesis. Since the 
log-area parameters are approximately ordered in terms of their sensi- 
tivity, the most expandable bits are those allotted to the lower-ordered 
g:’s.4 Depending on the speaker and the utterance, the bit rate can be 
lowered to between 900 and 1500 bits/s and still allow acceptable 
quality synthesis.* Figures 11, 12, and 13 illustrate the above examples 
for a bit-rate of 1400 bits/s (3; 3;1;4, 3, 2, 2,1, 1,1, 1,1,1,1). The 
synthetic speech in these examples is slightly degraded from the un- 
quantized synthesis, but the speech is still readily understood and the 
vocal attributes of the speaker are still apparent. It should be appre- 


* Acceptable quality speech synthesis is defined as speech containing all the in- 
formation content of the original without containing any annoying degradation in 
speech quality. 
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Fig. 8—2000-bits/s quantization of “Few thieves are never sent to the jug.” 
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Fig. 9—2000-bits/s quantization of ‘May we all learn a yellow lion roar.’’ (Female speaker, BG.) 
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Fig. 10—2000-bits/s quantization of “It’s time we rounded up that herd of Asian cattle.’’ (Male speaker, PB.) 
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Fig. 11—1400-bits/s quantization of ‘‘Few thieves are never sent to the jug.”’ (Male speaker, LG.) 
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Fig. 12—1400-bits/s quantization of “May we all learn a yellow lion 
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Fig. 13—1400-bits/s quantization of “It’s time we rounded up that herd of Asian cattle.’ (Male speaker, PB.) 
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ciated that the necessary storage requirements to produce acceptable 
quality synthetic speech in this method are nearly 3 the requirement 
for the pcm transmission of the LPc parameters (see Section I). 


V. REAL-TIME DPCM TRANSMISSION 


The prem scheme developed in Section III suffers from the draw- 
back that the calculation of the ppcm predictors and the quantizer 
step size are delayed until all the Lec parameters are available. For 
real-time speech synthesis, it is desirable that the process of parameter 
transmission be done concurrently with the measurement of the LPc 
parameters. In this section, we discuss several schemes for achieving 
real-time transmission while still retaining almost the performance of 
the optimum prem strategy. 


5.1 Average statistical system 


The first means of obtaining a real-time system is based upon the 
observation that the optimum pprcm first-order predictor for many of 
the Lpc parameters is nearly equal to one [b; = 1.0 in eq. (4) ]. Thus, 
the optimum linear prediction of the parameter x, is approximately 
x,-1. Table III is a comparison of the improvement factors Gop; ob- 
tained for b; = 1.0 and b, set equal to the optimum value. The overall 
improvement factors for b; = 1.0 are not significantly different from 
the optimum values, and the delay in calculating the optimum b, can 
be avoided by simply letting 6b; = 1.0. 

To design the optimum quantizer, it is necessary to know the stan- 
dard deviation of the signal to be quantized. However, our statistical 
studies have indicated that the standard deviation of the various 
difference signals are quite stable across different utterances and 
different speakers. Table IV shows the measured standard deviations 
for each difference signals computed with b, = 1.0. Table IV also 


Table III|— Comparison of G,,: in decibels with b, set equal 
to optimum value and b,=1.0. Sentence A is ‘Few 
thieves are never sent to the jug’ and sentence 
B is ‘May we all learn a yellow lion roar.” 


Pitch Power 91 92 
Sentence A 
b; = Optimum 23.7 12.2 14.7 12.2 
bi = 1.0 20.2 10.1 14.1 11.0 
Sentence B 
b; = Optimum 33.8 19.0 24.0 19.6 
b, = 1.0 33.1 18.8 23.9 19.2 
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Table !V— Measured standard deviations for the 
synthesis parameters 


Updated No Updating 

Pitch Period 13.01 16.5 

Power 27 X 105 27 X 105 
N 0.697 0.959 
g2 0.729 0.830 
g3 0.509 0.559 
94 0.510 0.554. 
gs 0.413 0.446 
96 0.417 0.430 
91 0.386 0.406 
gs 0.385 0.406 
Jo 0.377 0.399 
Ji0 0.346 0.364 
gu 0.332 0.342 
912 0.322 0.328 


contains the standard deviation for a system in which the prediction 
scheme is not updated for each unvoiced and voiced region. 

Using the standard deviations listed in Table IV and the quantizer 
discussed in Section IV, a robust transmission scheme is achieved. 
For example, the difference signal for the pitch period can be accurately 
quantized for differences as small as two samples or as large as 50 for 
three-bit quantization.* The synthetic speech quality for the average 
statistical system compares quite favorably to the optimum scheme, 
and has the added advantage of real-time implementation. 


5.2 Adaptive system 
5.2.1 Adaptive DPCM prediction 


The prom predictors can also be calculated without knowing the 
entire sequence of parameters by an adaptive method that is based 
upon the technique of “steepest descent.’’"! In this scheme, an initial 
estimate of the ppcm predictors is determined and then a new set of 
predictors is calculated to reduce the prediction error. The perturba- 
tion in the predictors is in a direction opposite the gradient of the pre- 
diction error taken with respect to the ppcm predictor vector. The 
resulting perturbation is given by 


5"(b;) = B-sgn (e,)-2,_;/ au |Ze—z |, (10) 
where B is the adaptation rate (typically, B = 0.09), and #, is the 


“Tf a nonlinear smoothing algorithm (Ref. 10) is applied to the raw pitch measure- 
ments, the variance of the corresponding difference signal is reduced by more than 
3. A two-bit quantization can then be used for pitch without diminishing the quality 
of the synthesis. 
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quantizer value of the parameter. For the prediction of the (r + 1)th 
sample of the parameter, the ppcm predictors are given by 


by +? = bj + 87(b;). (11) 


For a quantizer with B = 2, it can be shown that the adaptation 
scheme will match the changing spectral properties of the speech signal 
and result in near-optimum performance.” For the two methods given 
above, it should be noted that, in addition to the on-line calculation 
of the ppc predictors, it is unnecessary to transmit the predictors. 


5.2.2 Adaptive quantization 


In the previous section, the quantizer was constructed to take ad- 
vantage of the known properties or average statistical properties of 
each parameter’s difference signal. In this part of the paper, we intro- 
duce an alternate technique for estimating the signal variance. This 
method is based upan an adaptive approach developed by Cummiskey, 
Jayant, and Flanagan." In their scheme, a simple uniform quantiza- 
tion of the difference signal is used, but the step size for every new 
input is varied by a factor depending on which quantizer slot was 
occupied by the previous sample. Numbering the quantizer slots in 
the manner shown in Fig. 14, the updated step size A,41 is calculated 
from the previous step size A, by 


Arti = A,-M(|H,|), (12) 


where H, = 1, 2,---, B and the multiplier function M( ) is a time- 
invariant function of the quantizer slot number. 


OUTPUT 


INPUT 





Fig. 14—Numbering of quantizer slots for adaptive quantization. 
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To adequately match the step size to the signal variance, the mul- 
tiplier function must be properly chosen. Table V shows the multiplier 
functions found to be experimentally optimum for quantizing speech 
waveforms. Using this adaptive scheme (Appcm) and these multipliers, 
the quantization of the difference signals can also be efficiently achieved 
even when the initial step size is a poor estimate of signal variance. 
Table VI is a comparison of the signal-to-noise ratio for the adaptive 
scheme with a crude initial estimate of step size and the optimum 
quantizer discussed in Section IV. The results in Table VI are an 
encouraging demonstration that it is not necessary to know the sta- 
tistical structure of the difference signal to efficiently quantize the 
signal. In fact, it can be shown that, if the properties of the signal are 
nonstationary, the adaptive method is more suitable than the scheme 
used in Section IV. 

It should be noted that the above scheme does not apply for one-bit 
quantization (B = 1). A simple strategy for one-bit quantization has 
been developed by Jayant. Let c, and c,_; denote the values of suc- 
cessive bits in a one-bit scheme, then 


A, = AyiP rer, (13) 


where P has the typical value P = 1.5. Although this method was 
developed for quantizing speech waveforms, it performs quite well 
in quantizing the parameter difference signals. A comparison of this 
method and the optimum technique is shown in Table VII. Again, the 
adaptive scheme works well even with a poor initial estimate of signal 
variance. 


5.3 Synthesis 


To subjectively evaluate the performance of the adaptive methods 
suggested in this section, several speech utterances were synthesized. 
The synthesis scheme was again the one described by Atal and Han- 
auer,! but an adaptive quantizer and a second-order adaptive predic- 


Table V — Step size multipliers for B = 2, 3, and 4 (Ref. 7) 


2 3 4 

M(1) 0.80 0.90 0.90 
M(2) 1.60 0.90 0.90 
M(3) 1.25 0.90 
M(4) 1.75 0.90 
M(5) 1.20 
M(6) 1.60 
M (7) 2.00 
M(8) 2.40 
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Table VI — Comparison of the signal-to-noise ratio for the 
adaptive quantizer with crude initial estimate of step 
size and the optimum gaussian signal uniform 
quantizer. The analysis is for the sentence, 

“May we all learn a yellow lion roar.” 





gi G2 J3 
Bits | | 
Adaptive | Optimum | Adaptive | Optimum | Adaptive | Optimum 
2 12.6 13.6 18.3 18.4 15.6 16.5 
3 18.0 20.4 21.8 21.8 19.2 20.0 
4 22.8 21.9 24.9 23.9 24.0 23.1 


tion ppcm technique was used to transmit the Lpc parameters. The 
initial estimates of the predictors were b; = 1.0 and by = 0.0. A second- 
order analysis was performed because adaptive prediction makes the 
Gopt function saturate at a larger value than a nonadaptive predictor.’ 
The initial estimate of the quantizer step size was set equal to the 
standard deviations of the parameters listed in Table IV. For param- 
eters in which the quantizer uses only one bit, a first-order system 
with b; = 1.0 was used. 

Employing the bit assignment cited in Section IV, the quality of the 
synthetic speech was only slightly worse than the optimum scheme. 
Figure 15 shows a comparison of one example of the optimum scheme 
and the adaptive method. To achieve the performance of the optimum 
scheme, it has been found necessary to allot approximately one bit 
more per frame to the most sensitive parameters (usually pitch and 
power). 


VI. ORTHOGONAL LINEAR PREDICTION 


In the ppcm method of transmission, the value of the parameter 2, 
is predicted from previous values of the given parameter. However, 


Table VIl— Comparison of the signal-to-noise ratio for a 
One-bit adaptive quantizer and optimum one-bit 
gaussian signal uniform quantizer 


gi g2 Js 





Adaptive Optimum Adaptive Optimum Adaptive Optimum 


7.3 8.5 8.8 9.9 8.2 9.7 
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Fig. 15—Adaptive quantization (2000 bits/s) of “Few thieves are never sent to the jug.’’ (Male speaker, LG.) 
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the Lpc parameters have been experimentally determined to be quite 
redundant.!® Thus, the parameter x, can be predicted not only from 
its own past values but also in some sense from the values of the other 
Lec parameter. A more efficient method of transmission can then be 
obtained by exploiting all the available information about a given 
parameter. 

One means of exploiting the redundancy among the Lrc parameters 
is to generate a set of orthogonal parameters that are linear combina- 
tions of the original set. The new parameters are uniquely (one to one) 
related to the Lpc parameters and are calculated to be independent of 
each other and therefore do not contain any mutual information. If 
the original parameters are redundant, only a small subset of the 
orthogonal parameter will demonstrate any significant frame-to-frame 
variation. The process of obtaining the appropriate orthogonal param- 
eters is referred to as an eigenvector analysis.!® The orthogonal param- 
eters are termed eigenvectors, and each vector’s statistical variance 
is termed the eigenvalue of the eigenvector. 

To determine the eigenvectors, we first calculate the covariance 
matrix of the log-area parameters R across the utterance. If we denote 
gi; as the zth log-area parameter in the jth frame, then the elements 
of R are 


1 N 
rie = TD, (ges ~ MIs — ma), 
where 


1 X 
mn Bm 


and N is the number of frames in the utterance. Given the covariance 
matrix, the set of eigenvalues \,; are found by solving the set of simul- 
taneous equations 

[R — AI] = 0, 


where I is the identity matrix and |A| denotes the determinant of the 
matrix A. The eigenvectors ®; are then found as solutions of the 
equation 

A; — R@,. 


To illustrate the behavior of the Lpc parameters and the correspond- 
ing orthogonal parameters, Table VIII contains a listing of the typical 
variance (elgenvalues) of each calculated eigenvector parameter across 
the four utterances examined. The redundancy in the original log-area 
coefficients is reflected in the fact that more than 90 percent of the 
total statistical variance is contained in the first five or six eigenvectors. 
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Table VIIi— Measured eigenvalues for the four 
sentences analyzed: 
A. Few thieves are never sent to the jug. 
B. May we all learn a yellow lion roar. 
C. It’s time we rounded up that herd of Asian cattle. 
D. Should we chase those young outlaw cowboys? 





> 
ee) 
OQ 
o 


1 2.62 2.23 1.75 2.75 
2 1.44 0.80 1.29 1.58 
3 0.67 0.54 0.85 0.6 
4 0.44 0.38 0.52 0.36 
5 0.25 0.32 0.31 0.28 
6 0.21 0.24 0.17 0.22 
7 0.10 0.12 0.13 0.16 
8 0.09 0.10 0.10 0.15 
9 0.08 0.08 0.08 0.09 
10 0.06 0.05 0.06 0.06 
11 0.03 0.04 0.04 0.06 
12 0.02 0.01 0.02 0.03 


The higher numbered orthogonal parameters have a relatively small 
variance and can therefore be considered essentially constant through- 
out the utterance. Thus, the total information in the 12 log-area 
parameters can be effectively represented in the space of only the first 
six eigenvectors. 

The redundancy in the uPc parameters is not surprising in view of 
the fact that the speech signal can be synthesized with only three 
formant parameters (/";, F2, f;). Thus, the information contained in 
the 12 log-area coefficients are effectively duplicated in the space of 
only three formant parameters. The method of orthogonal linear pre- 
diction can be viewed as a constraint technique for squeezing the 
original parameters into a smaller but informationally equivalent set 
of parameters. The informationally equivalent set is formed by the 
most significant orthogonal parameters (significance is measured in 
terms of the standard deviation, or eigenvalue, of the orthogonal 
parameters). 

Experimental studies of a variety of speech utterances have shown 
that quite acceptable quality synthesis can be generated by trans- 
mitting only the six most significant orthogonal parameters, pitch, 
and power. The synthesis is performed by calculating the Lec param- 
eters from the transmitted orthogonal parameters and a prior: knowl- 
edge of the average values of the least significant orthogonal param- 
eters. For acceptable quality synthesis, only 22 bits/frame are needed. 
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The allotment of bits was as follows: 


Pitch: 3 bits/frame 

Power: 3 bits/frame 

Unvoiced-voice: 1 bit/frame 

First orthogonal parameter: 4 bits/frame 
Second orthogonal parameter: 3 bits/frame 
Third orthogonal parameter: 3 bits/frame 
Fourth orthogonal parameter: 2 bits/frame 
Fifth orthogonal parameter: 2 bits/frame 
Sixth orthogonal parameter: 1 bit/frame. 


The total transmission storage requirement in this technique is 
1100 bits/s for the synthesis parameters, 100 bits/s for the ppcm in- 
formation, and an initial one-time investment of 240 bits* for the 
necessary eigenvector information. Figures 16 to 18 illustrate the 
synthetic speech spectrograms generated by this technique for the 
examples previously examined. Depending on the speaker and the 
utterance, the bit rate for the synthesis parameters can be reduced to 
between 600 and 1000 bits/s and still yield acceptable quality speech. 
The low bit rate required for orthogonal linear prediction is quite 
attractive, but unfortunately this method involves a complex eigen- 
vector analysis and a delay in transmission to collect the statistical 
data necessary for the calculation of the eigenvectors. 


Vil. SUMMARY AND CONCLUSIONS 


The goal of this paper was the development of a more efficient 
method of transmitting the Lpc parameters. One proposed method 
involved the use of ppcm techniques. In ppcm transmission, we take 
advantage of the predictability of the parameter from its previous 
values to develop a more effective transmission scheme. Acceptable 
quality synthetic speech can be generated with prcm by allotting 
between 1000 and 1500 bits/s. This rate of information transmission 
is significantly better than the bit rates necessary for the conven- 
tional pcm methods. 

To enhance the practical application of the ppcm system, the 
methods of adaptive quantization and adaptive prediction were dis- 
cussed. These methods allow the on-line calculation of the ppcm 
predictors and quantizer step size. To further decrease the storage re- 


“Four bits for the average value of each of the six least significant parameters 
(24 bits) and three bits for each of the 12 coefficients required to compute each 
orthogonal parameter from the log-area coefficients (86 X 6 = 216 bits). 
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Fig. 16—Quantization of eigenvectors (1200 bits/s) of ‘Few thieves are never sent to the jug.” (Male speaker, LG.) 


1720 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1975 


YACOOOA NOILOIGSYd-YVANIT 


LeZL 


UNQUANTIZED 





QUANTIZED 








Fig. 17—Quantization of eigenvectors (1200 bits/s) of ‘“May we all learn a yellow lion roar.’”’ (Female speaker, BM.) 
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Fig. 18—Quantization of eigenvectors (1200 bits/s) of “It’s time we rounded up that herd of Asian cattle.’’ (Male speaker, PB.) 
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quirements of the Lec vocoder, the redundancy of the log-area param- 
eters was exploited. By transmitting only the most significant eigen- 
vectors, a considerable saving in bit rate can be achieved. 

The techniques discussed in this paper are not limited to the trans- 
mission of the Lpc parameters, but can also be used in conjunction 
with other vocoder systems. For example, the bit rate of a formant 
vocoder! can be reduced using a ppcm scheme for transmitting the 
necessary information. These transmission techniques have wide 
application and can prove very beneficial in a variety of synthesis 
schemes. 
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On the Theory of Self-Resonant Grids 


By |. ANDERSON 
(Manuscript received May 27, 1975) 


An approximate theory is developed to predict the frequency response 
of a self-resonant grid. The grid is comprised of capacitive and inductive 
elements and exhibits a band-stop resonance. The analysis is based upon 
the derivation, from physical considerations, of an equivalent circutt 
representation of the grid structure. Predicted results compare well with 
measured data. 


I. INTRODUCTION 


Arnaud and Pelow™ have recently described measurements of the 
transmission properties of several new types of self-resonant, metal grid 
structures. These grids, which are readily fabricated by photolitho- 
graphic techniques, have applications as millimeter-wave quasi-optical 
filters, or diplexers, in communications satellite antennas and in beam 
waveguide systems. The grid elements are symmetrical such that the 
grids may be used with two orthogonal polarizations. In this paper, 
we derive theoretical expressions for the frequency response of the 
simplest of the new grids and compare the results with measured data. 

The grid to be considered here is a periodic array of “‘Jerusalem”’ 
crosses as shown in Fig. la. We wish to determine the grid frequency 
response in terms of the dimensions of the elements when the planar 
transmitted wave is incident normally. On account of the complex ge- 
ometry of the grid elements, an exact treatment as a boundary value 
problem would be prohibitively difficult. Computer-oriented, numeri- 
cal techniques?? have provided a powerful means of solution for grid 
structures in the form of arrays of rectangular, or circular, apertures. 
The successful application of these techniques requires,* however, 
considerable caution in approximating the unknown aperture fields. 
When the aperture geometry is complicated, as here, this aspect of 
the numerical approach poses a considerable difficulty. 

Now, in general, the transmission properties of grid structures can 
be described® in terms of an equivalent impedance, together with a 


“In eq. (2) of Ref. 1, \ should be replaced by /2. 
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(a) 


Fig. 1—Jerusalem-cross array and approximate equivalent circuit. 


section of transmission line which represents propagation in free space. 
For example,® consider the transmission of a plane wave incident nor- 
mally upon a grid of thin, perfectly conducting, parallel metal strips of 
period p. When p <j, where 2 is the wavelength, the equivalent 
impedance is a shunt inductance, or capacitance, depending upon 
whether the electric vector of the incident wave is parallel to, or per- 
pendicular to, the edges of the strips. In the following section, an 
approximate circuit representation of the present grid is derived from 
physical considerations and from the known results for grids of 
parallel strips. This approach lends itself to a simple understanding of 
the grid transmission properties and, furthermore, leads to useful 
design formulae. 


Il. ANALYSIS 


As shown in Fig. la, the period of the array is p, the width of the 
inductive strips is w, and the separation between adjacent crosses is 
g. The length and width of the capacitive segments of each cross are 
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d and h, respectively, and the thickness of the grid is ¢. It is assumed 
that 
t<K<w<p, hKp <i, and gKd<K d., (1) 


The electric field, #‘, is incident normally on the grid with the electric 
vector directed as shown. For purposes of discussion, we shall refer 
to this as the ‘‘vertical’” direction; the incident magnetic field, H’, is 
then in the horizontal direction. The effect of the vertical “‘dipoles,”’ 
each of length d and width A, at the sides of the crosses is negligible 
for d < X. It is therefore assumed that only current that flows parallel 
to EH‘, along the vertical inductive strips and across the horizontal 
capacitive strips, is significant in determining the field scattered by 
the grid. On the basis of this assumption, we now consider the mag- 
netic and electric fields in the vicinity of the grid. 

Since w <p and h <p, the magnetic field about the grid, due to 
current flowing along the vertical inductive strips, is approximately 
the same as that about a corresponding uniform inductive grid of 
period p and strip width w. Hence, the stored magnetic energy of the 
Jerusalem-cross grid may be represented approximately by the equiva- 
lent inductive reactance, X(w), of this uniform grid, where’ 


X(w) = c {i  eosee (=) + FQ, ¥)| (2) 
and é 
7 Ce pe 2 
PO, u) = 7 P2oa+| Ra - sf, (3) 
with 


a=[1-(2) "1; ¢ = cost (52), s=1-—c. (4) 


The reactance X(w) is normalized with respect to the intrinsic im- 
pedance of free space. The first term in (2) can be derived’ from mag- 
netostatic considerations; the second term is a correction factor which 
is negligible when p < X. Since t < w, the effect of thickness upon the 
inductive reactance is negligible. 

With regard to the distribution of electric field, it is noted, from 
symmetry considerations, that there is no component of electric field 
normal to the grid on the planes A and A’ of Fig. la. Without dis- 
turbing the electric field we may, therefore, insert a pair of infinitely 
thin, perfectly conducting plates at A and A’ which are perpendicular 
to the plane of the grid and distance p apart. In the quasi-static case, 
when p<), the electric flux about the grid elements within this 
parallel-plate transmission line is concentrated between the gaps of 


SELF-RESONANT GRIDS = 1727 


the horizontal capacitive segments. We assume this concentration to 
be maintained at all frequencies for which p < 2. Since g <d, the 
effect of fringing at the extremities of the segments is negligible and the 
electric flux, per unit width of the parallel-plate line, is d/p times that 
of a corresponding uniform capacitive grid of period p, gap width g, 
and thickness ¢. This implies that the stored electric energy, of the 
Jerusalem-cross grid, may be represented approximately by an equiva- 
lent capacitive susceptance 


Bg, 1) = < Bulg t), (5) 


where B,(g, ¢) is the (normalized) susceptance of the corresponding 
uniform grid. For the case t = 0 we have® 


B.(g, 90) = #2 tin | cosee (#2 ) | + F(A, of, (6) 


where F(\, g) is given by (3) with w replaced by g. The equivalent 
impedance of a uniform capacitive grid of thickness ¢ includes® a seg- 
ment of transmission line of length ¢. When ¢ < 0.5), this transmission 
line may be represented by a IJ-network of shunt capacitors and a 
series inductor. In the present case, t <\, the series element may be 
neglected and the total susceptance is!* 


B,(g, 1) = Bug, 0) + ae. 


(7) 


The second term in (7) may be derived equivalently by considering 
the additional (parallel-plate) capacitance introduced by the finite 
thickness of a capacitive diaphragm in a parallel-plate transmission 
line of height p. From (5), (6), and (7), the capacitive susceptance of 
the Jerusalem-cross grid is approximately 


BG) = “ jin | cosee (#2) | 4+ F(A, 9) + x (8) 


We have obtained approximate values of reactances with which to 
describe the stored magnetic and electric energies of the grid and now 
consider the equivalent circuit representation. It has been assumed 
that only current that flows vertically along the inductive strips, and 
across the gaps of the horizontal capacitive segments, is significant 
in determining the transmission properties of the grid. This suggests 
that the Jerusalem-cross grid can be represented approximately by a 


“In Ref. 10, the sign of the second term for B,; in eq. (83) on p. 200 should be 
positive. 
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reactance, X,, where 


1 

Ag AM) Bi aayS (9) 
shunted across a transmission line of impedance Z, as shown in Fig. 1b. 
The impedances in the equivalent circuit are normalized with respect 
to the impedance (Z,) of free space, and the impedance seen by a plane 
wave incident normally on the grid is Z;. The grid transmission co- 
efficient is now expressible in terms of the grid reactance X,. The input 
impedance, Z,, is 

jXo 


T+ 9X, oY 
The corresponding voltage reflection coefficient, R, is 

Zi-1 
sare 

and the grid power transmission, | 7'|?, is 

4X5 
a, (ree Prima oe 

IT? = 1 — [R[ = > P (12) 


Substituting (2) and (8) into (9) and (12) then gives the grid trans- 
mission response in terms of its dimensions. To the present order of 
approximation, |7'|? is seen to be independent of h when h <p. 
A first-order approximation for the rejection wavelength d,, defined 
by the equation X, = 0, is readily found by assuming p <) so that 
the terms F(A, w) in (2) and F(A, g) in (8) may be neglected. Further- 
more, if the effect of grid thickness is also neglected, by putting ¢ = 0, 
and if the cosecants are replaced by the small argument forms, we find 


po P 2p 

X(w) 2 In (22), p< (13) 
~, 4d 2p _ 

Bg, 1) & > in (22), p<r, t= 0. (14) 


From (9), the wavelength, \,, at the rejection resonance is then 


coe 2p 2p \. 
he 2 ap in (22) in (22) (15) 


This result provides an approximate functional dependence of the 
resonant wavelength upon the grid geometry. 

The effect of the dielectric sheet which supports the grid has been 
neglected in the preceding analysis. In general, the presence of an 
adjacent, low-loss, dielectric layer will increase the grid susceptance, 
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Fig. 2—Transmission response of Jerusalem-cross array. 


B(g, t), by modifying the electric field in the vicinity of the capacitive 
gaps. In the case examined by Arnaud and Pelow,! however, the sheet, 
which has a relative permittivity of about 2.5, is thinner than the grid 
itself and, as such, is not expected to modify the grid transmission to 
the present order of approximation, 


It. COMPARISON WITH MEASUREMENTS 


The grid measurements of Arnaud and Pelow! were conducted under 
approximately plane wave conditions and for a range of incidence 
angles from 5 to 45 degrees.* It was found that the frequency of the 
rejection resonance, and the shape of the transmission response, were 
practically independent of the angle of incidence within this range. 
Figure 2 shows the predicted frequency response for normal incidence, 
as obtained from (12), compared with measured data for an incidence 
angle of 5 degrees. The experimental curve is from Fig. 3 of Arnaud and 
Pelow’s paper and is for a grid of dimensions p = 1.400 mm, d = 0.750 
mm, w= h=0.180 mm, g = 0.090 mm, and ¢ = 0.018 mm. The 
shape of the transmission response is predicted well by the theory; 
the error in the prediction of the rejection frequency is 7 percent. The 
first-order expression (15) for the rejection wavelength is within 10 
percent of the value obtained from (12). 


IV. CONCLUSIONS 


We have examined the transmission properties of a self-resonant 
grid that is comprised of capacitive and inductive elements. An 


“No measurements were taken at exactly normal incidence to avoid multiple 
reflections within the measuring system. 
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approximate theory has been developed to predict the frequency 
response of the grid when illuminated by a plane wave at normal inci- 
dence. The theory is based upon the construction of an appropriate 
equivalent circuit in which the values of the reactances are obtained 
by modification of known solutions for simple, parallel strip grids. A 
comparison of results with measured data shows an error of 7 percent 
in the prediction of the grid rejection resonance. By way of compari- 
son, the corresponding approximate expressions (2) and (6), for parallel 
strip grids as obtained from rigorous analyses,’ can be in error by 
about 1 to 5 percent over the range of frequencies considered here. 
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An Approximate Method for Calculating 
Delays for a Family of 
Cyclic-Type Queues 


By S. HALFIN 
(Manuscript received May 22, 1975) 


A study of the marker-register dial-tone delay problem in No. & crossbar 
switching machines led to a special type of cyclic queuing model. In this 
paper, we present a method for calculating approximately the steady-state 
delays of an arriving customer. When applied to the marker-register 
problem, the model emphasizes the order in which markers are assigned 
to waiting calls and the fact that part of the markers’ time ts unproductive 
when an “all registers busy’’ condition occurs. Some numerical results are 
presented, which agree with the observed phenomenon that, for a constant 
marker load, the delays of Touch-Tone® calls are influenced by the load 
on the dral-pulse originating registers, and vice versa. The results are 
compared to those of a simulation of the same problem. The numerical 
results compare favorably in the range of loads that produce a dial-tone 
speed of between 0.05 and 0.18. 


I. INTRODUCTION 


The queuing model described in this paper resulted from a study of 
the marker-register dial-tone delay problem in No. 5 crossbar switching 
machines. A number of queues with Poisson arrivals of equal rates 
are served in a cyclic order by a server with constant service time. 
Upon arriving at a nonempty queue, the server chooses a customer 
from the queue at random. After one service time, the customer either 
leaves the system with a certain predetermined probability or rejoins 
his queue. In both cases, the server uses a fixed amount of time and 
moves to the next queue; thus, at most one customer leaves the system 
following each arrival] of the server at a queue. 

Related models were treated by Cooper,! Cooper and Murray,” and 
Eisenberg.® In Refs. 1 and 2, the server either empties the queue being 
served or serves all those present at the queue in its arrival epoch. 
The case of two queues with different arrival rates is treated in Ref. 3. 
In these papers, the Laplace-Stieltjes transforms of the waiting time 
distributions were obtained. Attempts to obtain the distributions for 
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similar models by approximate methods were made by Leibowitz‘ and 
Schay, Jr.® 

The method presented here approximately calculates the steady- 
state delays of an arriving customer. The approximation is carried out 
by modifying the model to achieve a manageable state space. Next, 
the method is applied to the problem of dial-tone delay in the No. 5 
crossbar switching machine. The order in which markers are assigned 
to waiting calls and the fact that part of the markers’ time is unpro- 
ductive when an “‘all registers busy” condition occurs are emphasized 
in the model. The numerical results presented agree with the observed 
phenomenon that, for a constant marker load, the delays of Touch- 
Tone® calls are influenced by the load on the dial-pulse originating 
registers, and vice versa. The results are compared to those of a simu- 
lation of the same problem. 

Several features of the No. 5 crossbar machine, which may have an 
influence on the dial-tone speed, were excluded. Some of these features 
were investigated in subsequent work by H. A. Guess® and are de- 
scribed in more detail in Section XVI of this paper. 


ll. THE MODEL 


Let N queues Fy, Bo, ---, Hy, N = 2, be given. Customer arrivals 
to each queue constitute independent Poisson processes, all having the 
same rate \. The queues are served by a single server in the following 
way: At each point in time, the server is at some queue. Transitions 
in its position occur at discrete time epochs which are equally spaced 
with periods of duration 7. At such time epochs, the server moves 
instantaneously to the next nonempty queue in a circular order, chooses 
a waiting customer in that queue at random, and stays with this 
customer until the end of the period. The served customer then leaves 
the queue with probability p*, or remains in the queue with probability 
q* = 1 — p*. If all the queues are empty at a transition point, then 
no change occurs in the position of the server. 

We assume further that 7’ is small with respect to the accuracy 
with which we want to know the delays, and thus all the arrivals can 
be assumed to occur at the transition epochs and the queuing process 
can be considered in discrete time. Thus, if X;,, is the number of 
arrivals to ; at the kth time epoch, then all X;,,,1 575 N,k 2 0, 
are independent identically distributed random variables, all having 
a Poisson distribution with mean AT. 


Ili, THE FULL STATE SPACE 


The state of the queuing system at any time epoch is defined as the 
(N + 1)tuple (m1, me, ---, my, n), where m; is the number of waiting 
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customers at /;, and n is the position of the server before the transition. 
We call the set of all such states the full state space. It is then clear 
from the discussion in the previous section that our queuing system, 
with the full state space, is a stationary Markov chain. Thus, in 
principle, one can calculate transient and steady-state probabilities. 
However, the full state space is much too large for practical computa- 
tional purposes. 

Consider, for instance, the dial-tone delay problem where a typical 
number of N is 15. Then, even if the system is so underloaded that 
each mm; can be restricted to be either 0 or 1, we have 2)° & 15 ~ 500,000 
states. A natural approach, which we follow in the remainder of the 
paper, is to approximate the behavior of our system with systems 
having a smaller size state space. 


IV. THE “BLACK BOX” APPROACH 


Some important facts about the system can be deduced by consider- 
ing only the total number of customers in all the queues s = my, + 
m2. + --- + my. It is clear that the system with this single state is 
again a stationary Markov chain. It is, in fact, a discrete version of 
an M/D/1 queue, with the added feature that a customer who is held 
by the server returns to the queue with probability q*. Alternatively, 
the service time measured in units of 7 may be considered as having 
a geometrical distribution with mean 1/p*. The traffic intensity of 
the system is the p = ANT'/p*; thus, we have Theorem 1. 


Theorem 1: A necessary and sufficient condition for the nonsaturation 
of the system ts NNT < p*. 


Note that, because of the symmetry, a particular queue in the sys- 
tem is saturated if and only if the system as a whole is saturated; 


hence, Theorem I provides a saturation condition for all the individual 
queues. 


Let us denote by A; the probability that a Poisson-distributed ran- 
dom variable with mean \NT attains the value 7. The system has the 
following transition probabilities. For s > 0, 


Pr (s > 8’) = Ayesarp* + Ay—sq* 
and 
Pr (0— 8’) = Agr. 


The equations for the steady-state probabilities P, are then 
s’+1 
Py = AyPo+ DY) (Ag-spip* + Av-g)P. 8 = 0,1,---. (1) 
s=1 
Equations (1) can be solved recursively, starting with Py) = 1 — p. 
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One can also calculate the generating function 


>, Pus 
s=0 

(pe) (a) 
ug* + p* _ yer -u) ) 


P(u) 





where a = ANT. By evaluating P’(1), the expected value of the total 
number of customers in the system § is 


2a — a? 


gaa) (2) 
V. THE MODIFIED MODEL 


We now consider our original system with a new state space. A 
state will now consist of a triplet (m, M, n), where m = m1, M = me 
+ms3 +--+ + my, and n is the same as previously, namely the posi- 
tion of the server before a transition. It is clear that the new state 
space is much smaller than the full state space; however, the Markovian 
property is lost. This can be seen by the following argument: If I 
was positive at time & — 1, and if in the transition between k — 1 tok 
the server skipped a large number of queues, then those M customers 
were concentrated in the remaining queues; thus, it is probable that, 
in the k to k + 1 transition, a small number of queues will be skipped. 

At this point, we modify our model to make it a Markov chain with 
respect to the new state space. To do this, we need to define one-step 
transition probabilities so that the behavior of the modified model 
will approximate that of the original model. Let the position of a 
customer be the queue number where he waits. Our key assumption 
concerns the probability distribution of the positions of the M cus- 
tomers, given the state (m, M, n). 

For the remainder of this section, let us enumerate the queues 
Ij, +--+, Hy by starting with the queue following the position of the 
server and observing the cyclic order, skipping /;. Next, we make the 
following assumptions: Let M > 0; then 


(2) The positions of the M customers are independent, identically 
distributed, random variables. 

(it) The probability that any one of the M customers will be in the 
7th queue (in the new order) is 7:;(/), where 


(N —7i—1)R(M) 
M(N — 2) 


where 6(M) = 1/(N — 1) — [R(M)/2M] and determination 
of R(M) is described below. 


wi(M) = b(M) + 





24,9. tN ad, 
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The rationale behind assumption (77) is that a customer is less likely 
to be in a queue that has just recently been visited. The average differ- 
ence between the number of customers in the first and last queues is, 
according to assumption (72), M[mi(M) — ay_i1()] = R(W), which 
should be approximately equal to the expected number of arrivals 
during one full cycle of the server. Thus, 


R(M) ~ XTR,(M), 


where Ro(J) is the expected number of nonempty queues, given VM. 
Finally, we approximate Ro(M) by 


where the right-hand side is the expected number of nonempty queues 
among E., ---, Hy, if the M customers are uniformly distributed. 
We conclude this section by using the new assumptions to calculate 
some probabilities that will be needed later. 
Let J denote the number of successive empty queues following the 
position of the server (when £ is disregarded), J = 0,1, ---, N — 1. 
The distribution of J depends on M. Let 


Q;(M) = Pr 2 7 — 1) 
qi(M) = Pr (J = 3 — I) 
j 1, 2, ed Vs 


Using assumptions (z) and (22), we have 


N-1 


Q:(M) = ee r(t |" $5 bo ea 


0 if M>O0 
Qv(M) = 
1 if M=0, 


and 
qi(M) = Q;(M) — Qj41(M) jg=1,---,;N-1 
qn(M) = Qn(M). 


VI. TRANSITION PROBABILITIES FOR THE MODIFIED SYSTEM 


Given a state (m, M, n), the transition probabilities of the position 
of the server can be expressed in terms of the Q;(J/) and q;(M)’s as 
follows. 
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For m > 0, 
dnt—n(M) if n’>n 
Pr(n-> 7’) = 4Qneni(M) if n’=1 
0 iL <n Sn. 
For m = 0, M > 0, 
Qn-n(M) if n’>n 
Pr (n— n’) = 40 if 7’ =] (3) 
QN+n'’—n-1 if 1 < n! < Nn, 
and, form = 0, M = 0, 
Dae aS an 
Pr (n> 7’) = 
0 if n’#n. 
Let us denote by A; and a; the probability that Poisson-distributed 
random variables with means \(N — 1)T and XT, respectively, attain 
the value 7. We can now write the state transition probabilities. 


Form >0,M>0: 


Pr (im, M, n) > (in’, M’, n’)] 
Ant —mQnt'—n(M) CA aa ip* + Aweug*] if n’>n 
= Aa —mMQnva1.—n(M) dint myrp* + fneng” | if ni = i 
0 otherwise. 


For m = 0, M > 0: 


Pr [(0, M, n) — (m', M’, n’)] = Qw[Ase—acyip* + Aar—ug*] 
dnt—n(M) if n’>n 
x 40 if n’=1 


dn+nt—n—1(M) if 1 < n’ nN. 


IA 


For m > 0, M = 0: 


Pr [(m, 0, n) — (m’, M’, n')] 
Aw’ LQm'—m41p* —+ Gen ge" | if n’ = l, 


0 We se 1, 
and finally, form = M = 0, 
Om Au U n’=n 
Pr [(0, 0, n) > (m’, M’, n')] = (4) 
0 if n' ~n. 


Vil. STEADY-STATE EQUATIONS FOR THE MODIFIED SYSTEM 


If we consider the total number of customers in the system s = m 
+M, then it is clear that all the results of Section IV remain valid 
for the modified system. In particular, ANT’ < p* is the necessary 
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and sufficient condition for the existence of a steady state. In the 
following discussion, we assume that this condition is satisfied. Let 
P(m, M, n) be the steady-state probability of the state (m, M, n) 
satisfying the following equations. 

For n’ > 1, m’ 2 0, M’ = 0: 


Pim’, M’, n’) 


n’—1 om’ M’4+1 
= om Am'—m|Yn'—n (M)[A aer—aryip* + Aw—uq* | 
n=l m=0 M=1 
N M’41 
x Pn, M, n) + Ds 2» Am’ A mr—ar4ip* nie Aar—_uq* ] 
n=n’ M=1 


Xdnin'—n—1(M)P (0, M, n) + Om Aw P (0, 0, n'). (5) 
For n’ = 1, m’ 2 0, M’ = 0: 


M* 


N m’+4+1 
Pim’, M’,1) = x yz Amr—mQwyy1—n(M) [ane mysp* a Gnicnd | 
n=l m= 0 


XP(m, M,n) + am AP (0, 0, 1). 


M= 


Vill. NUMERICAL SOLUTION OF THE STEADY-STATE EQUATIONS 


Equations (5) can be solved either as a system of linear equations 
with the auxiliary equation 
> P(m, M,n) =1 
m,M,n 
or by starting with any initial distribution and iterating it through (4) 
until a desired degree of convergence is obtained. 

It seems adequate to adopt the second method since the steady state 
of the modified system is of interest to us only as an approximation 
to the steady state of the original system. Thus, an extensive compu- 
tational effort to obtain an accurate solution to (4) is not warranted. 
A good initial distribution to start the iterations can be obtained in 
the following way. We have 


N 
Pn Mon) =P, s=0,1,-- 


m+M=s n=1 


*) 


where the P, were computed by the method outlined in Section IV. 
If we divide the s customers uniformly among the queues and make 
the position of the server random, we get the following distribution : 


1 1 
Po(m, s — m,n) = wn P (m, =, s) P, 
s=0,1,:-- m=0,-°',8 n=1,---,N, 
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where B(m, 1/N, s) is the binomial probability of m successes out of 
s trials, with probability of success = 1/N. 


IX. EVALUATION OF THE DELAYS FOR THE MODIFIED SYSTEM 


Let f(m, M,n,k) be the probability that a customer arriving at 
Ei, will wait k service periods before leaving the system, given that 
on his arrival the system went into the state (m, M, n). The customer 
stays for at least one service period, sok = 1. Notice also that m = 1. 


Theorem 2: The delay probabilities satisfy the following recursive 
equations: 


* 
f(m, M,n, 1) = © Qyanyx(M) 





and 
fim, M, n,k +- 1) = 0, + 02 
m2z1, M20, l1sn&SQN; k21, 
where 
om Omt—mQnt—n(M) [A aerate rp* + Am _—muq* | 
n'=n+l M’=M-—1 m!/=m 
op ad Xf(m', M', n’, k) if M>0,n<N 
0 if M=0O0orn=WN 
02 = = x Am—mQnii-n(V) ae ip* + Am'—mg* 
ie ps — +l—n m m!—m-- m'—m 


xX f(m'’, M’, 1, k). 


Proof: f(a, M,n, 1) equals the probability that the server moves to 
£;, that the particular customer is selected for service, and that he 
leaves the system after the service period. Hence, the formula for 
k= 1. 
For k = 1, we have 
f(m,M,n,k+1)= D> Prim, M, n) > (m’, M’, n’) 
m’,M’',n’ 


(\ the customer stays in £y]f(m’, M’, n’, k). 


Using eq. (4), we get that o, is the part of the right-hand side cor- 
responding to n’ # 1 and oz is the part of the right-hand side corre- 
sponding to n’ = 1. 

Theorem 2 provides a method for calculating the delays conditional 
on the state. Let f*(k) be the probability that a customer arriving at 
fi; will wait k service periods before leaving the system, given that 
before his arrival the system was in the steady state. 
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Theorem 8: 


C.) ce) N 
Fe= LV DV LY fim’, Mn’, k)P*(m', Mn’) kb = 1,2, °°, 


m'=1 M'=0 n/=1 


where 


—_ 





n’'—1 M'+1 
P*(m', M', nn’) = 


P(m', M',n')-— a0 DD anr—-n(M) 
1 — ao n=1 M=1 


(Aareeargip* + Aar—oug*)P(m’', M, n) | 


and 


1 N M’' 
P*(m’, M’,1) = ee | Pon’ M’,1) — ao x Pa Qya1—n(M) 


[p*P(m! + 1, M, n) + g#P(m’, M, n)3| 
form’ 21,M’20,2Sn' SN. 


Proof: The theorem is valid if it is shown that P*(m’, M’,n’) is the 
probability of the state (m’, M’,n’) at the point of arrival of the 
customer, say, point u, for all m’ 21, M’20, 1S n'’ SN. The 
probabilities of states with m’ = 0 is zero, since there is at least one 
customer in £;. We know that at u — 1 the system was in a steady 
state. The transition probabilities between u — 1 and wu, conditional 
that at least one customer arrives at £1, are obtained from (4) by re- 
placing a; with aj, where 
1 


* 1 — ao 
ai= 





a; i=1,2,::: 


0 7= 0. 


The expressions for P*(im’, M1’, n’) can now be calculated by operating 
on the steady-state probabilities with the modified transition prob- 
abilities and using the fact that the steady-state probabilities satisfy 
eqs. (5). 


X. VALIDITY OF THE MODIFIED SYSTEM 


The difference between the original and the modified systems is in 
the rules of the server movement. In the modified system, the server 
does not follow the cyclic order. However, to calculate the delay dis- 
tribution, we are interested only in the pattern of the time points 
when the server is in £;, and that, hopefully, is similar to the cor- 
responding pattern in the original system. The degree of similarity 
is difficult to check, except by simulating the original system. Verifica- 
tion of the “‘reasonability’’ of the modified system may be made by 
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checking whether each FE, receives the same number of visits by the 
server, and if /, gets the right amount of expected number of waiting 
customers, 1.e., 


3 Pon wie n=1,-::,N 
m=0 M=0 N 


and 


0 00 N - 
Lm dv YY P(m,M,n~s, 
m=0 N 


M=0 n=1 


where § is given by eq. (2). 


XI. THE MARKER-REGISTER SYSTEM 


In the rest of the paper, we apply the model to the dial-tone delay 
problem in the No. 5 crossbar switching machine. Following are some 
operational features of the No. 5 crossbar switching machines which 
are relevant to the dial-tone delay distribution. 

Calls appear on line link frames (LiF). The dial tone markers (ptm) 
which are not busy are paired to the waiting calls. Under ‘“normal”’ 
operation, i.e., when several pTMs are free, the LuFs look for available 
DTMs according to a fixed preference order. When all the prMs become 
busy, a gate is closed and the pTMs serve first those LLFs that contain 
waiting calls at that moment. If an LLF has more than one call waiting, 
only one call will be served during the gating period. 

When a pTMm becomes idle following the ‘‘all markers busy’”’ condi- 
tion, it looks for a waiting call according to the following scheme. Each 
pTM has its own order in which it scans the uirs. Thus, for example, 
when there are four ptMs and 60 LLFs, the first ptm will scan the LLFs 
in the natural order from 0 to 59, the second ptm will start at Lur 15, 
go to 59, and then come back to 0 to 14, etc." 

When a pr locates an LLF with waiting calls, it chooses one of those 
calls and proceeds to look for an originating register (or) for the call. 
The above choice may be considered random for all practical purposes.' 
If the ptm finds a vacant or, then it connects the call to the or, and 
the calling customer gets a dial tone. If no or is available, then the 
pTM releases the call, and it continues to wait in its LLF and to bid for 
a pTM. In both cases, the holding time of the pTm is constant and 
approximately equal. We denote this time by 7’, where T is approxi- 
mately 0.25 second. In fact, this time is approximately 0.21 second in 


* This is the recommended arrangement, although not all No. 5 crossbar entities 
observe it. 


T We omit consideration of the systematic preference for serving calls in vertical 
group 2. 
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the case where no or can be found, but we ignore this difference to 
simplify our model. 

We also assume that the distribution of holding times of the ors 
is negative exponential for a conservative estimate of the delay dis- 
tribution. The arrival of calls to each LLF is assumed to be Poisson, 
with the rate being equal for all LLFs. We denote the rate for a single 
LLF by X. Finally, in many cases there are two types of calls, dial-pulse 
and Touch-Tone, where both types are served by the same pTMs but 
require different ors. When there are two types of calls, the ratio be- 
tween their arrival rates is assumed to be the same in every LLF. 


XII. A QUEUING MODEL FOR ONE TYPE OF CALL 


The system described in the previous section is quite complicated, 
and it appears that, to model such a system and be able to derive 
numerical results from the model, some simplifying assumptions are 
inevitable. One such model was proposed by W. 8S. Hayward.’ Its basic 
assumption is that, in order to be served, a call must find both a 
marker and a register idle. Once the marker and register start process- 
ing a call, they act independently of each other, each having exponenti- 
ally distributed holding times. To solve the resulting state equations, 
Hayward introduced a system with one type of server, which approxi- 
mates the behavior of his model. 

The present queuing model emphasizes the order in which the 
markers are assigned to waiting calls, and takes into account the 
fact that the time a marker spends serving a call is nearly the same, 
whether or not it found a free register. 

First we assume that each pTM serves only those LLFs which are of 
high priority on its list. Thus, in the example of the previous section, 
the first ptm will serve only the first 15 LL¥s, the second pTM will serve 
only the next 15 LuFs, etc. Such an assumption is Justified under heavy 
load conditions. We denote by N the number of LLFs which are served 
by one DTM. 

Next, we assume that each pT serves its LLFs in a cyclic order and 
that, whenever it finds a LLF with waiting calls, it serves exactly one 
call. This assumption is asymptotically valid under heavy traffic 
loads, because of the gating procedure described in the previous 
section. 

Finally, we assume that whenever a DTM serves a call there is a 
fixed probability p* that an or will be available and thus that the 
waiting time of the call will end (i.e., the customer gets a dial tone). 
This assumption would hold if the availability of the ors is independent 
of the number of waiting calls, which is clearly not the case. This 
assumption will cause our model to somewhat underestimate the 
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delays, while the first assumption tends to overestimate them. The 
value of p* can be taken approximately as the delay probability of a 
call, given that the arrival process of the calls to the ors is Poisson, and 
therefore it can be computed by the Erlang C formula. 

Thus, we arrive at the model which was described in Section II, 
with the server the ptm. The server chooses a customer (call) from 
the queue at random. After one service time, the customer either leaves 
the system with probability p* or rejoins his queue (that is, waits in 
his LLF). The server then moves to the next LuF having a waiting call. 


XIII. SOME THEORETICAL AND NUMERICAL RESULTS 


It was shown in Section IV that the occupancy of the pr is 


5 = AND. 
p* 





Hence, a necessary and sufficient condition for nonsaturation of the 
system is ANT < p*. Also, the expected total number of waiting calls 
in the LLFs which are served by the prM is 


§= a ay (2 — PP”): 


Thus, for a fixed occupancy of the prm, the expected total number of 
waiting calls is a monotone decreasing function of p*. The same is 


true for the expected delay, W, since by Little’s formula 


— s $ 
WS ag ape 


= 1 2 
W = saa): 

A standard measure for the quality of service is the dial-tone speed 
(pts), which is the probability that the call will have to wait three 
seconds or more for a dial tone. Figure 1 presents some computed 
values of the prs for various values of p* and \ with p held constant 
at three different values. It is seen that prs is also a monotone de- 
creasing function of p*, for a fixed pTM occupancy. 





and so 


XIV. A MODEL FOR TWO TYPES OF CALLS 


Consider now the case of a system having dial-pulse and Touch-Tone 
calls; this is the usual situation in No. 5 crossbar offices today. Let 
the arrival rates from each LLF be \1 and de, and let the probabilities 
of finding available registers be p; and 73 for dial-pulse and Touch-Tone 
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0.64 NO. OF DTMs = 4 
NO. OF LLFs = 60 


ise DTM HOLDING TIME = 0.25 SECOND 


0.48 


0.40 


0.32 


DIAL—TONE SPEED 


0.24 


0.16 


0.08 


p= 0.84 
X= 0.224p* 





0.25 0.50 0.75 1. 
p* (AVAILABILITY OF REGISTERS) 


Fig. 1—pts as a function of p* and A for constant DTM occupancy p (queuing model). 


calls, respectively. The loads on the ptm due to the two types of calls 
are 


_ MNT 
eee | 








_ MNT 


and pe 
Po 


To approximate the delays of this system, we consider a system with 
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one type of call, for which \ = Ay + Az and p = pi + p2. The appro- 
priate p* of this system satisfies 


Piss A 8, 

— Ma/pi + 2 /p2 

Now the state probabilities can be computed by the approximate 
method; however, it is necessary to modify the formulas for the 
computation of the delay distributions to obtain those distributions 
conditional on the type of call. 

Let fi(m, M, n, k) be the probability that a call of type z (¢ = 1, 2) 
arriving at LLF No. 1 will wait & service periods before leaving the 
system, given that on its arrival the system went into the state 
(m, M, n). The recursive formulas of Theorem 2 have to be modified to 


Pp 


fo(m, Mm, 1) = Tt Qy—nea(M) 


and 
fi(m, M, n, k - 1) = 061 + oe 
ie, MS Iepen:. 224 


where o; is as in Section IX except that f is replaced by f; and 


m — 
m 





1 * 
Am'—m+1)P 


eo ie} 
a= yi » 
m'=M m'=m— 


(Aatr-aiQat-n(M) 


We di Ue aps 
+ Om'—m ( m q a m ) [osc ’ M ’ 1, k) ]. 
The proof of the validity of the modified formulas is along the same 
lines as the proof of Theorem 2. 


XV. NUMERICAL RESULTS 


Several computer runs were made for a typical large system with 
60 LLFs, 4 ptMs (N = 15), 100 dial-pulse ors and 50 Touch-Tone ors, 
with both dial-pulse ors and Touch-Tone ors having a mean holding 
time of 13 seconds. 7’ was taken to be 0.25 second in all runs. The 
parameters varied were \ = \i + 22, the total input rate per LLF, and 
a = \i/de (the ratio of the rates of the two types of calls). 

Figures 2 and 3 describe the results for a = 2, that is, when the 
ratio of the loads is the same as that of the ors. Figure 2 describes the 
behavior of the occupancies of the ors and the pro, while in Fig. 3 
the prs is plotted as a function of \. Figures 4 and 5 present the cor- 
responding results for a = 3. The values of \ were chosen to be near 
the point of saturation, i.e., where the occupancy of the DTM ap- 
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Fig. 2—Occupancies of the ptm and ors as a function of \ (queuing model). 


proaches 1. In this range, the pts is sensitive to small perturbations 
in A. 

Figure 6 shows the dependence of the pts on a for a fixed \ (A = 
0.158 per second). It can be seen that the quality of service deteriorates 
as a diverges from the neighborhood of the ratio of the number of 
dial-pulse ors/number of Touch-Tone ors (which equals 2 in our case). 
This is consistent with the observation that Touch-Tone delays are 
significantly influenced by the dial-pulse or load for a constant offered 
load to the pTM. 

The results were compared with those of a simulation model which 
we constructed for the system described in Section XI. Tables I and 
II present a comparison between the results of the queuing model and 
the simulation. We compare the intensities of input for which levels 
of the pts are reached between 0.05 and 0.25 for a balanced system 
and between 0.05 and 0.15 for an unbalanced system. Examining those 
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Fig. 83—prts as functions of \ (queuing model). 


Table |— Values of A for which given levels of DTS are reached, 
for Ai/ dg = 2 
(a) Dial-Pulse Calls 
DTS Queuing Model Simulation 
0.05 0.1680 0.1695 
0.10 0.1705 0.1735 
0.15 0.1720 0.1760 
0.20 0.1730 0.1775 
0.25 0.1735 0.1790 


(b) Louch-Tone Calls 


DTS Queuing Model Simulation 
0.10 0.1690 0.1690 
0.15 0.1705 0.1720 
0.20 0.1720 0.1740 
0.25 0.1725 0.1760 
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Fig. 4—Occupancies of the ptm and ors as a function of \ (queuing model). 


Table Il — Values of A for which given levels of DTS are reached, 


DTS 


0.05 
0.10 
0.15 
0.20 


DTS 


0.05 
0.10 
0.15 


for Ay/A2 = 3 


(a) Dial-Pulse Calls 


Queuing Model Simulation 
0.1540 0.1550 
0.1565 0.1590 
0.1580 0.1610 
0.1585 0.1620 

(b) Touch-Tone Calls 

Queuing Model Simulation 
0.1575 0.1600 
0.1590 0.1625 
0.1600 0.1640 


CYCLIC 
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Fig. 5—prts as functions of A (queuing model). 


tables, we observe that the differences in the corresponding figures 
for the two models in these regions are less than 4 percent of the total 
input. Also, it can be observed that the computed pts grows faster 
in the queuing model than in the simulation. The reason is that the 
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Fig. 6—pts and ptm occupancy as functions of Ai/A2 for a fixed A, A = 0.158. 


assumption made in Section XII, that the probability of finding all 
registers busy can be computed by the Erlang C formula, is incorrect 
in this region. Figure 7 compares the pts as computed by the two 
models. Again we conclude that the fit is fair in the “critical’’ region. 
It would have been desirable to validate the results by performing a 
field trial. Such a trial] should consist of measuring the pts for a No. 5 
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Fig. 7—pts as function of \ (queuing model and simulation). 


crossbar machine with high constant loads. However, observing Fig. 8, 
where the cumulative distributions of the hourly pts from the simu- 
lation were plotted, one sees that these distributions have long tails. 
This implies that, to get a reasonably accurate estimate for the pts, 
say, with a standard error of 1 percent, we would have to run the 
trials for around 25 hours while keeping the load constant. This seems 
to be a difficult task. 


XVI. SUMMARY AND CONCLUSIONS 


We presented a method for modifying the original model, as pre- 
sented in Section II, to a model which has a much smaller state space. 
Methods were described for calculating the steady-state distributions 
of the states and of the delays in the modified model. The model was 
applied to the problem of calculating dial-tone delays in the No. 5 
crossbar switching machine. This was accomplished by making some 
simplifying assumptions about the order of service of the waiting calls 
by the markers. The numerical results were compared to those of a 
simulation, and found to be close on an important range of the pts. 
This gives us a certain amount of confidence that both models are 
valid, which is especially important because of the difficulty in validat- 
ing the models by experimental data, as discussed in Section XV. 
However, the reader should be aware that several features of the No. 
5 crossbar machine, which may have an influence on the DTS, were 
excluded. 
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Fig. 8—Cumulative distributions of the hourly dial-tone speed (simulation). 


A major effect not covered by the present model is the effect of 
horizontal group blocking on dial-tone speed and on dial-tone marker 
waste usage. Recently obtained field data and theoretical studies re- 
ported in a subsequent paper by H. A. Guess® have shown that dial- 
tone speed and dial-tone marker occupancy can be appreciably in- 
creased by horizontal group blocking caused by high average line link 
frame loads and also by poor load balance. Consequently, the dial-tone 
speeds associated with a given call origination rate in an actual No. 
5 crossbar office may be higher than would be predicted by our model. 
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Recently obtained field data and theoretical studies show that, for a fixed 
subscriber calling rate, dial-tone delay in No. 5 crossbar can be appre- 
ciably increased by high average line-link frame loads and also by poor 
load balance. The increased delay is caused by waste dial-tone-marker 
usage generated by a small number of calls that encounter horizontal 
group blocking in obtaining a dialing connection. This paper discusses 
an analytical model to predict the time from receiver off-hook to receipt of 
dial tone under various service conditions in No. 6 crossbar. 


l. INTRODUCTION 
1.1 Description of the dial-tone connection process 


In a No. 5 crossbar switching machine, dial tone is provided to a 
subscriber line, terminating on a line-link frame (LLF), by an originat- 
ing register (oR), terminating on a trunk-link frame (TLF), via a series 
of three network links: line links, junctors, and trunk links. The dialing 
connections are set up by dial-tone markers (pTMs), which are common 
control devices. Each line-link frame contains a number of crossbar 
switches that are used to establish connections between subscriber 
lines and trunks, or between subscriber lines and service circuits, such 
as originating registers. The crossbar switches that form line concen- 
trators on which groups of subscriber lines terminate are called hori- 
zontal groups. Maximum size offices typically contain from 4 to 6 
ptMs, 40 to 60 LLFs, 20 to 30 TLFs, and up to 140 ors. Each LiF contains 
10 horizontal groups and each horizontal group is a concentrator con- 
taining between 29 and 59 subscriber lines on the input side of the 
switch and 10 line links on the output side of the switch. 

To provide dial tone to a subscriber, the off-hook signal from the 
subscriber line initiates a bid for a ptm through a connector circuit 
unique to each uur. As soon as a DTM becomes available, it locates an 
unoccupied or and then attempts to find a dialing path (consisting of 
a line link, a jJunctor, and a trunk link) connecting the or with the 
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subscriber line. As soon as the connection is established, the prm 
releases and proceeds to serve other calls waiting for dial tone. The 
OR provides dial tone, receives the dialed digits, obtains a completing 
marker, transmits information to it, and releases. The completing 
marker then establishes the connection between the calling subscriber 
and an outgoing or intraoffice trunk. 

If all ors are busy, the ptm releases and the call rejoins what is 
effectively a queue of calls waiting for a ptm. If an or is available, but 
no (unoccupied) dialing path connecting the or with the subscriber 
line can be found, the pt will release the or, obtain a second (usually 
different) or, and try to find a dialing path between that or and the 
subscriber line. When a ptm is unable to find a dialing path between 
a given or and a given subscriber line, a matching failure is said to 
occur. If, on the second try, the ptm cannot find a dialing path, the 
DTM releases and the call rejoins the queue of calls waiting for a pTm. 
In such a case, 3, DTM second-failure-to-match (DTM2FTM) is said to 
occur. 

The method of assigning DTMs to waiting calls is controlled by a type 
of “gating” circuitry designed to equalize service and to reduce the 
incidence of long delays in obtaining dial tone. When one or more 
pTMs are free (light traffic operation), the LLFs look for DTMs according 
to a fixed preference order. When all ptms become busy and a request 
for a DTM occurs, a gate is closed and only the LuFs containing calls 
waiting for dial tone at that moment are put inside the gate. The 
pTMs then proceed to serve the LLFs inside the gate. If more than one 
call is waiting on an LLF, only one call will be served during the gating 
period. Once a call on an LLF is served during a given gating period, 
that LLF is put out of the gate, whether or not the ptm is successful 
in establishing a dialing connection for the call. When all LuFs with 
requests at the start of the gating cycle have been put out of the gate, 
the gate opens; if there are sufficient waiting calls to cause all of the 
pTMs to become busy again, a new gating cycle will be started; other- 
wise, light traffic operation will resume. During a gating period, a ptm 
that becomes idle scans the uuFs in cyclic order. Each ptm has a 
different starting LLF for the scan so as to equalize service. (The de- 
scription of the gating procedure is taken from Refs. 1 and 2.) 


1.2 Effects of matching failures on dial-tone delay 


Repetitive matching failures can occur in establishing the dialing 
connection essentially because the first-stage crossbar switch on which 
the subscriber line terminates (the horizontal group) is a concentrator 
whose output links (line links) have holding times that are much 
longer than the holding times of the common control devices that set 
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up the dialing connections (the prms). The average line-link holding 
times, being largely determined by conversation holding times (with 
allowance for ineffective attempts), can be on the order of 150 seconds 
or more, while the ptm holding times are typically on the order of 
0.25—-0.40 second. 

Although matching failures (blocking) can occur when one or more 
of the 10 line links on a horizontal group are unoccupied, repetitive 
matching failures under such conditions are rare because a DTM is 
quite likely to be successful in establishing a dialing connection after 
a few attempts. The expected holding time in a blocked condition of 
a call that finds all 10 line links busy is the lesser of the length of time 
for one of the 10 line links to become free (with a short added time 
for DTM uses and matching failures in setting up the dialing path) and 
the length of time that a subscriber is willing to wait. 

Since, at average busy-hour load levels, less than about one percent 
of all originating calls are predicted to encounter an all-10-line-links- 
busy condition, previous dial-tone-delay studies have assumed the 
effect of matching failures on dial-tone delay to be small. However, 
the fact that line-link holding times are much longer than pra holding 
times means that a call which finds all 10 line links busy can remain 
blocked for long enough to consume a large number of DTM uses (except 
at calling rates sufficiently high that very little of the offered blocked- 
call load is carried by the ptms). Thus, it is possible for a small number 
of calls experiencing blocking to generate a disproportionate number of 
waste DTM uses, increase DTM occupancy, and thereby increase dial- 
tone delay for all other calls in the office. 


Il. ANALYTICAL MODEL (LIMITING FORM) 
2.1 Assumptions 


Since the main effect of matching failures on dial-tone delay is 
caused by the resulting waste pTmM use, and since the dial-tone delay 
distribution of the small proportion of calls experiencing repetitive 
matching failures can be calculated approximately (see Section 5.2), 
the effect of matching failures on the dial-tone delay distribution of 
calls that do not experience repetitive matching failures has been repre- 
sented in terms of a queuing model with two classes of calls, good calls 
and bad calls, defined as follows: 


(2) A good call experiences no matching failures but is subject to 
delay caused by pTMs and ORs. 

(2z) A bad call experiences total network blocking (no dialing con- 
nection available) and defects from the system after an expo- 
nential waiting time. 
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We will first describe the mathematical structure of the analytical 
model used for dial-tone delay calculations. We will refer to this model 
as the limiting model. Next, we derive the prm saturation load* and 
prove that it is not changed by the presence of bad calls. 

In Appendix B, we prove that the equilibrium queue length and wait- 
ing time distributions of the limiting model are the limits in distribu- 
tion of a sequence of equilibrium distributions arising from a model in 
which the expected bad-call arrival rate approaches zero and the 
expected bad-call waiting time (until defection) approaches infinity 
in such a way that their product, total erlangs of bad calls, is constant. 

In the limiting model, a good call arrives and finds a random 
(truncated) Poisson-distributed (but time-independent) number of 
bad calls permanently present in the system.' The queue discipline 
is characterized by random order of service. The prms cannot dis- 
tinguish between good and bad calls when choosing a call to be served. 
This corresponds to the fact that, in the No. 5 crossbar switching 
machine, a DTM cannot recognize that a 10-line-links-busy condition 
exists on a particular horizontal group and, hence, cannot avoid wasting 
time serving calls for which no dialing path exists. Equilibrium good- 
call queue length distributions are computed conditional upon the 
number of bad calls present in the system. Calculating the expectation 
of the conditional distribution over the distribution of bad calls gives 
the unconditional] dial-tone-delay distribution for good calls. The con- 
ditional distributions depend on the total office calling rate and on the 
number of LLFs, DTMs, and ors but do not depend on the horizontal 
group load or the load variation. Hence, delay distributions for a range 
of frame load and balance effects can be calculated using the same set 
of conditional distributions, thereby greatly reducing the computer 
time needed for parametric studies. 

Since ptm holding times are approximately constant (for a given 
set of office parameters and traffic characteristics) and since these 
holding times are small with respect to the accuracy with which it is 
necessary to be able to predict delays, we treat the dial-tone delay 
process as a discrete time queue with a constant service time of T 
seconds. Good calls are assumed to arrive in batches according to a 
Poisson process at times k7, for k = 1, 2,---. Immediately upon 
arrival, the good calls join the queue of good and bad calls waiting 
for dial tone. Calls are chosen at random from the queue for service 


“The ptm saturation load is defined to be the good-call originating load below 
which a steady state good-call-queue-length distribution exists and above which 
such a steady state distribution does not exist. 

T When the bad-call input is assumed to be peaked, the number of bad calls in 
oe vaien has a (truncated) negative binomial distribution. This is discussed in 

ection IV. 
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by the ptms with each call—whether a good call or a bad call—having 
an equal probability of being chosen. A good call served on the ptm 
cycle beginning at time kT will either acquire an or, receive dial tone, 
and thereupon exit the system at time (k + 1)7, or else will fail to 
acquire an or and will return to the queue of calls waiting for dial 
tone at time (kK + 1)T. A bad call served on the ptm cycle beginning 
at time kT will return to the queue of calls waiting for dial tone at 
time (k + 1)T. Note that since new arrivals occur only at the time 
points kT, a pr that is idle on the cycle beginning at time k7’ will 
remain idle at least until the start of the cycle beginning at time 
(k + 1)T. The dial-tone delay for a call is the length of time from the 
moment when the call arrives to when it obtains dial tone and (simul- 
taneously) leaves the dial-tone queue. Thus, the minimum possible 
dial-tone delay in the model is 7 seconds. 

In the analytical model, two prms serving N line-link frames are 
used to represent four pTMs serving 2N line-link frames. A single pT 
holding time equal to the office average ptm holding time is used for 
both good and bad calls. In actuality and in the simulation models, 
bad calls have somewhat longer pr holding times than good calls 
and good calls that encounter a condition of all-ors-busy have some- 
what shorter ptm holding times than good calls that do not encounter 
a condition of all-ors-busy. Comparison of results from the analytical 
model with those of the gating simulation model indicates that these 
simplifications tend to offset each other in the range of interest. 

A further simplification in the analytical model concerns the manner 
in which availability of ors is treated. Since No. 5 crossbar offices 
typically contain over 100 ors and since the ors are in tandem with 
the pras, it is presently not practical to keep track of the number of 
occupied ors directly in an analytical model. Availability of ors is 
treated by assuming that at each time point kT’, all ors are busy with 
probability gq and two or more ors are free with probability p = 1 — q. 
The calculated probability that exactly one or is free and that a 
dialing connection is available between this one or and the given sub- 
scriber line is sufficiently small, in the occupancy range of interest, 
so as not to warrant the additional complexity caused by introducing 
this effect. 

The idea of using a discrete time model and the method of treating 
or availability through use of a fixed probability of all ors busy are 
due to Halfin.* Following Halfin, the probability q is taken to be the 
erlang C’ probability of all ors busy. At or occupancies above about 
0.90 with frame loads low enough that few second failures to match 
occur, this method of treating or availability somewhat underpredicts 
the delay caused by an all-registers-busy condition (based on com- 
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parison with simulation results). The underprediction arises because, 
at higher or occupancies, once all ors become busy they tend to remain 
busy for a time period equal to several ptm cycles, as one would expect 
from the erlang C delay formula. 

Preliminary studies were made using an analytical single-server 
cyclic queuing model, developed by S. Halfin,*? which represented the 
No. 5 crossbar gating process in considerable detail but did not take 
into account the effects of horizontal group blocking. These studies 
showed that the delays predicted by the cyclic queuing model do not 
differ appreciably from those of a discrete-time M/D/1 queue with 
feedback and random order of service. The latter model requires 
about 1/30th of the computer time required by the former. Both 
models overpredicted simulations. For these reasons, explicit repre- 
sentation of the gating process was not attempted and a discrete- 
time M/D/2 queue with feedback and random order of service was 
taken as the starting point for developing an analytical dial-tone-delay 
model to include effects of horizontal group blocking. 


_-— RANDOM ORDER OF SERVICE 
SIMULATION MODEL 


aoe ANALYTICAL MODEL 


oc GATING SIMULATION MODEL 


PROBABILITY THAT DIAL—TONE DELAY EXCEEDS X SECONDS 





SECONDS OF DIAL—TONE DELAY 


Fig. 1—Comparison of analytical and simulation models. Dial-tone delays are at 
1400 ccs/uuF and 0.70 ptm occ (excluding pTM2FTMs). 
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Comparison of simulation results with predictions using the limiting 
analytical model (based on two markers and random order of service) 
shows that the predicted delays typically fall somewhat above or very 
close to simulated delays based on four markers with gating order of 
service and somewhat below simulated delays based on four markers 
with random order of service. In the light of these results and because 
of the large scatter in actual measured No. 5 crossbar dial-tone delays, 
it did not seem worth the considerable added complexity to include 
explicit representation of the gating process in the analytical model. 
Typical results are shown in Figs. 1 and 2. 


2.2 Queue length equations for the limiting model 


The queue-length process for good calls in the limiting model is a 
discrete-time Markov chain with finitely many irreducible classes 
{C,,} that are noncommunicating in the sense that no transition be- 
tween different classes is possible. Hence, each class is itself an irre- 
ducible Markov chain. The kth class is the queue-length Markov chain 
for a discrete-time M/D/2 queue with random order of service, 


_— RANDOM ORDER OF SERVICE 
x SIMULATION MODEL 


--GATING SIMULATION MODEL 


= 


~~ ANALYTICAL MODEL 


PROBABILITY THAT DIAL-TONE DELAY EXCEEDS X SECONDS 





SECONDS OF DIAL—TONE DELAY 


Fig. 2—Dial-tone delays are at 1600 ccs/LLF and 0.80 ptm occ (excluding DTM2FTMs). 
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“feedback” (occurring when an all-ors-busy condition is encountered), 
and k blocked calls residing permanently in the system. The pTMs 
cannot distinguish these & permanently blocked bad calls from the 
good calls in the system, so each of the k bad calls and each of the good 
calls in the system compete on an equal basis for pts. Upon comple- 
tion of service by a pT, a bad call rejoins the queue. With probability 
q, the good calls served at any given time period rejoin the queue and 
with probability p they leave it. We will show that each of these 
queuing systems has the same saturation load and we will describe 
how the equilibrium-queue-length distribution and the equilibrium- 
waiting-time distribution for each system is calculated. 

Viewing the whole process again from the standpoint of a Markov 
chain with finitely many noncommunicating classes, the queue-length 
and waiting-time distributions referred to above may be regarded as 
being conditional on k, the number of blocked calls (permanently) 
present in the system. Now let k be a random variable with a truncated 
Poisson distribution of mean x, or, equivalently, regard the Markov 
chain as having any initial distribution 7(k, 7), where 7 denotes the 
number of good calls in the system and where the marginal distribu- 
tion «(k, -) of the number of bad calls in the system is truncated 
Poisson with mean x. Then the (unconditional) equilibrium-queue- 
length and equilibrium-waiting-time distributions for this system may 
be computed by taking the expectation (with respect to the truncated 
Poisson distribution of k) of the conditional queue-length and waiting- 
time distributions for each of the individual systems represented by the 
classes C’,.. 

Let X, denote the number of good calls in the queue at time n7' and 
let Y denote the number of bad calls (permanently) present in the 
queuing system. Let 


P,(i, 7) = Pr(Xeu = j|X, =tand Y =k) 


and let P,(¢) be the equilibrium-queue-length distribution for the case 
of k bad calls. Let A(n) be the probability that n good calls arrive in 
one service time interval (of length 7). Then, by assumption, 





A(n) = * fi 201 as 
=0 otherwise, 
where 
A= ANT 
\i = LLF originating call rate (calls per second on one LLF), 
N = number of tuFs served by the two DTMs. 
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The transition functions P;(-, -) are given by 
P(0, 7) = ACJ) 
Pod, 3) = pA(j) + GAG — 1) 
Pot, J) = pA(j —t +2) + qA(j7—7) fori 2 2, (1) 
and for k = 1, 
P..(0, 7) = ACJ) 
P(t, Jj) 
= pLk(k—1A (j—1)+20kA (7-741) +7(¢-1)A(G—-7+2) // 
[A+a)e+7-1)]+qAgG-a for 21. (2) 


Note that P,(7, 7) = Oforz > 7 + 2 and that P,(j + 2, 7) > 0. 
If the equilibrium distributions P;,(-) exist, then they satisfy the 
equation 


Pra) = & Puld)Pali, i). (3) 


Let f, be the generating function for the kth queue length process; 
then 


file) = = P, (ides. 


The generating function fo is easily seen to be given by 


pe plz — IE (2 + 1Po(0) + 2Po(1)]_ 


fo(z) zed-z) — (p + gz?) (4) 
If Po(-) exists, then lim.+, fo(z) = 1, and hence 
r 

Po(0) + 3Po1) = 1 - ap (5) 


Let S{(p) denote the saturation load of the kth process with two 
DTMs, as discussed above, and let S{?(p) denote the saturation load 
of the analogous one-pTm system with k blocked calls. It can be shown 
that S{?(p) = p for all k = 0, 1, 2, ---, ie., a necessary and sufficient 
condition that a steady-state queue length distribution exists in the 
one-DTM system is that \ < p.* Using this result, we will show that 
SY (p) = 2p for all k = 0,1, ---. 

Since the Markov chain for the kth process above is irreducible, 
P,(0) > 0 and P,(1) > O whenever the stationary distribution Po(-) 
exists. Hence, it follows from eq. (5) that S§?(p) S 2p. Since SP (p) 


* This result can be proved using a theorem of Kushner (Ref. 4) and some properties 
of generating functions. 
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< SE’ (p) S SP (p) S 2p for ky = k, = 0, it suffices to show that 
S(p) = 2p for all k = 0. 

Consider a two-pTmM system with 2k blocked calls, which differs from 
the 2kth process defined above only in that (7) each ptm maintains a 
separate queue with half the originating traffic being assigned to one 
of the ptmMs and the other half being assigned to the other pTm and 
(ii) each prm serves k blocked calls. Let So.(p) denote the saturation 
load of this modified system. Then clearly, S2.(p) S S#(p) and, by 
symmetry, So.(p) = 28{(p) = 2p, because the modified system 
simply consists of two identical one-marker systems working sep- 
arately, with each being assigned half the incoming traffic. Hence, 
S$? (p) = 2p and so S!(p) = 2p for all k = 0. 

Thus, for \/2p < 1, a stationary distribution exists for the kth 
process for each k = O. In the sequel, we will confine ourselves to the 
case where \/2p < 1. In this case, fo exists and is analytic for |z| <1 
and continuous for |z| = 1, so the numerator of the expression on the 
right-hand side of eq. (4) must vanish for any |z| S 1 for which the 
denominator vanishes. It is easy to see that the denominator has a 
single real root on the open interval (—1, 0) and that z = 1 is a root. 
By applying Rouche’s theorem to the functions 2? exp [A(1 — z) ] 
and p + gz along the circle |z| = 1 + «, one can show that for \/2p 
< 1 and for e > 0 sufficiently small (depending on \ and p), the ex- 
pression 2? exp [A(1 — z)] — (p + qe?) = 0 has exactly two roots in 
the open dise |z| <1-+ ¢«. Thus, the root on (—1, 0) and the root 
at 1 are the only roots within the closed unit disc and 


(& + 1)Po(0) + EPo(1) = 0, (6) 
where £ is the unique root of the expression 
aeri-2) — (p + gz”) = 0, ot We Oe ae <a 8 (7) 


Solving eqs. (5) and (6) for Po(0) and Po(1) yields the unique 
solution 


(1 — £)p 
_ (1+ &)(2p — 2) 
se a (l—é)p ’ 8) 


where £ is defined by (7). The values Po(7) for 7 2 2 can now be com- 
puted recursively from (8) using (8). 

When k > 0, computation of P;,(-) by the method of generating 
functions involves numerical solution of some rather unwieldy 
differential equations. The following approach, which makes use 
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of properties of recurrent Markov chains, is easier to implement 
computationally. 

In Appendix A we show that, for \/2p <1, Px(-) is given by 
c(a + 76), where c and r are constants and a, 6 are the two (particular) 
eigenvectors, defined by (9) below, of the transition matrix P;(2, 7).* 
We further show that the condition, a; + r8; 2 0 for each j, deter- 
mines 7 uniquely and provides an easy way to calculate r. Once r is 
obtained, the constant c is uniquely determined by the condition that 
> j20 Px(7) = 1. Since the quantities a; and 8; can be computed re- 
cursively, the above observations lead to a simple algorithm for com- 
puting P,(-). The results are given below, the proof is in Appendix A. 
Let 


apo = ] 
ay = 
n+1 
an — 25 ajP,(j, n) 
aay RCRD od ea > 
An+2 Pita n=O (9a) 
Bo = 0 
Bi = 1 
n+1 : 
Bn — 2 BP uC, ) 
ee oe > 
Bn4e Pi(n as 2, n) n = 0. (9b) 
It is shown in Appendix A that there is a unique constant r satisfying 
ee ae | alln = 3 (10a) 
lim m, =r = lim M,, (10b) 


where, for n = 3, the increasing sequence m, and the decreasing se- 
quence J/, are defined by 
Na = max | — Be <0, 6; > 0| 
j 


jsn 


(11) 


I 


M, rnin | — “Fa; > 0, B; <0}. 


jin B; 
Since the quantities m, and J7, can be computed recursively, eqs. 
(10a), (10b), and (11) yield a well-defined algorithm for computing r. 
Once r has been computed, P;(7) can be calculated from (3), setting 
P,(0) = c, Px(1) = re and determining c by the requirement that 
Lj=0 Px (7) = 1. 


* Each of ¢, r, a, and @ depends on k. For each integer 7 2 0, a; and 6;, denote the 
jth components of the eigenvectors « and 8, respectively. 


DIAL-TONE DELAY 1765 


Since recursive calculation of P;(7) by (3) involves successive sub- 
tractions, it is essential to perform all computations in double precision. 

In a No. 5 crossbar switching machine, the measured DTM occupancy 
during an hour is defined to be the total work time (in seconds) of all 
the ptms divided by the product of the number of ptms and the number 
of seconds in an hour. Thus, the measured occupancy is the average 
fraction of time that a given pT is occupied. In the mode] the (steady- 
state) DTM occupancy is defined to be the limit, as n approaches in- 
finity, of the expected fraction of time that a given pTm is occupied 
during the n work cycles of length 7 beginning at times 7, 2T, ---, nT. 
This limit is equal to the probability that a given ptm is occupied on 
a work cycle beginning with the system in steady state. Hence, 


DTM occupancy = P (two or more calls are waiting for dial tone) 
+ 4P (exactly one call is waiting for dial tone) 
= 1 — P (no calls are waiting for dial tone) 
— iP (exactly one call is waiting for dial tone). 


As discussed earlier, the number of bad (blocked) calls waiting for 
dial tone at any given time has a truncated Poisson distribution in the 
model. Let 


x = total erlangs of blocked calls on two prMs 
k = maximum possible number of blocked calls that can be in the 
system (waiting for dial tone) at any time* 
K yk 
CcK >= 2X i! . 


Then the probability that k blocked calls are in the system is given 
by (x*/k cre. Let p, denote the (steady state) prm occupancy for 
the case of no blocked calls in the system and let p» denote the (steady 
state) pTM occupancy including the effect of blocked calls. It follows 
from (5) that 

pu = d/2p. (12a) 
Hence, 
= 1 — czxPo(0) — 3 [erePo(1) + xereP1(0)] 
= 1 — cxx(1 — pu) — (@/2)ezRP1(0). (12b) 


zo) 
o 
| 


2.3 Delay equations for the limiting model 


In this section, we calculate the dial-tone delay probabilities for a 
call that arrives when the system (just prior to the arrival of the call) 


“The truncation parameter K determines the number of conditional delay dis- 
tributions to be calculated, so K should be taken to be no larger than needed to retain 
sufficient accuracy in the calculations. For values of x in the range of interest e~*¢zx 
differs from 1 by less than about 10-3 for K = 5; thus, for computational purposes, 
K may be taken to be 5. 
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is in steady state. Recall that the queuing model is in discrete time 
with arrivals at times n7’ and DTMs operating at times n7’. Calls served 
on the prm cycle beginning at n7' either exit the system at time 
(n + 1)T or else return to the queue at time (n + 1)T. 

Thus, the steady-state probability that a call experiences a dial- 
tone delay of (m + 1)T seconds can be computed recursively in terms 
of the probability of a delay of m7 seconds by straightforward condi- 
tional probability calculations involving matrix multiplications. 

Let 


Y = number of bad calls (permanently) present in the system 
X, = number of good calls in the queue at time nT 
X = number of good calls in the queue immediately after an arrival 
when, just prior to the arrival, the system was in steady state. 


Thus, X is the total queue length just after arrival of a call when the 
system is in steady state. Let 


Queue is in steady state at time (n — 1)T 
with k bad calls permanently present in| _ 


Pr(j) = P Xn = the system, and at least one good call 
arrives at time nT 
Then 
P[X = j/¥ =k] = Py (j). 
Let 


Pie ap | Xara = re =i, Y = k, and at least one or | 


call arrives at time nT 


and 


Aw) = 295 nz=1 


= 0 otherwise. 
Then A(n) is the conditional probability that n calls arrive at time 


nT given that at least one call arrives at time n7' and P,(i, 7) is defined 
by eqs. (1) and (2) with A used in place of A. Also, 


~ . gr? v\ PB ° * 
PCJ) = 2 PiOPiG, J). 
Let W.(z, k) be the conditional probability that a call arriving in 
steady state has a dial-tone delay of nT seconds, given that the queue 
length of good calls upon arrival of the call is ¢ and that the number 


of bad calls permanently present in the system is k. Thus, letting D 
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denote the dial-tone delay, we have 
W.G@,k) = P[D=nT|X=iand Y=k] 721. 
Then 


ee. 


: F 2 
W(t, k) = min 1, Ea) p 
and Wi4i(-,k) can be computed recursively from W,(-, k) as de- 
scribed below. 
Let X, = 7 2 1 and consider any one given call out of the 7 calls 
present at time nT’. Let 


Xnii = jand the given call does 
@,.(t, 7) = P| not leave the system on the ptm 
cycle beginning at time nT 





Xn =7 and y=} 
Then 
Worlds b) = Y Oxi, Wa (Gh), 
and, for z = 1, @:(2, 7) is given by* 
ov(i, f) = [1 —™ 29 |tpagj e+) + aA — 91 
+ [| EE? oa 0 


t 








oi = | 3 | BL@- MaG — 6+) 424-741 
+A - 9] +[ 35 Joau-0 
and, fork = 2, 
t—1 k 
ei) =[1- p55] ]? ae i Ne aaa 


Ce 
+A] + (Gog) au-a. 


The second term on the right-hand side of the above equation is 
the probability of the event that: (2) the given call zs selected for 


“The derivation of @x(i, 7) for k = 2 is given below. The derivations of @;(i, j) 
for k = 0 and k = 1 are analogous and are omitted. 


1768 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1975 


service by the pTMs, (72) all ors are busy, and (277) Xn41 = 7. (Recall 
that, in the model, either all ors are busy, with probability g, or else 
two or more ors are available, with probability p = 1— q). The 
first term is the probability that: (2) the given call is not selected for 
service and (72) Xnz41 = j. In the first term, the expression in large 
brackets is the conditional probability that X,+1 = j given that the 
call is not selected for service. This conditional probability is itself 
composed of the probabilities of two mutually exclusive events. The 
second term within the large brackets is the (conditional) probability 
of the event that all ors are busy and X,4; = j. The first term within 
the large brackets is the (conditional) probability that two or more 
ors are free and X,41 = 7. The sum from / = 0 to / = 2 in this term 
pertains to the cases where 0, 1, or 2 good calls, respectively, are selected 
for service by the pTMs. Since this sum is part of a probability condi- 
tional upon a given good call not having been selected by the prMs, 
the available population from which the prms may select calls con- 
sists of k bad calls (k = 2) and z — 1 good calls (¢ = 1). Since the 
selection is without replacement, the selection probabilities have the 
hypergeometric form shown above. 
By the law of total probabilities, 


PLD = nP|¥ =k] = > W.G, MPLS = 1¥ = 4), 


P[D > nT|Y =k] =1— X P(D=mTlY =k], 


m=1 


and 


K Cyne 
PLD > nT] = > SEF PLD > aT |¥ = ki], (13) 





where 


K 
Cox = d. 2*/kl. 
k=0 


Ill, CALCULATION OF OFFERED BLOCKED-CALL LOAD 


In the limiting analytical model, the number of blocked calls in the 
system has a time-independent truncated Poisson distribution with 
mean x. This section describes a method for computing z, the mean 
offered erlangs of bad calls, in terms of the distribution of carried load 
among horizontal groups in the office. Using the limiting analytical 
model in conjunction with these methods for calculating the offered 
blocked-call load, we can calculate the No. 5 crossbar dial-tone delay 
distribution in terms of the calling rate and the distribution of carried 
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load among horizontal groups in the office. Thus, we can predict the 
effect of frame load and balance on dial-tone delay in No. 5 crossbar 
offices. 

We first construct a model for an individual horizontal group and 
express the expected bad-call load from one horizontal group as a 
function of the carried horizontal group load. To get the total ex- 
pected bad-call offered load for the office, the expected contribution 
from an individual horizontal group is integrated (numerically) over 
the office horizontal group load distribution. The model, described 
below, for representing an individual horizontal group may be called 
a modified finite source Palm delay model. 

We assume that the input to a horizontal group is from a finite 
number of sources with equal calling rates and that call-holding times 
are exponential. The number of sources N is taken to be 35 although 
in actuality most horizontal groups have 49 or 59 subscriber lines. 
The reason for the use of the lower number of sources is that an earlier 
study by W. S. Hayward, Jr.® showed that blocking on concentrators 
with unequal line occupancies can be approximated by blocking calcu- 
lations based on equal calling rates and a lower number of sources. 
The use of 35 sources was suggested by J. G. Kappel.® 

The calculations take into account the fact that an incoming call 
cannot occupy a subscriber line when all ten line links are occupied. 
Calls that find all line links busy will either defect or will eventually 
obtain a line link. While waiting for a line link to become available, 
a call is assumed to have an exponential waiting time until defection, 
with a mean of 380 s.* 

In the case of ideal load balance, each horizontal group in the office 
is assumed to have a true carried load of Zz erlangs. In the case of less 
than ideal load balance, the distribution of true carried load among the 
individual horizontal groups is assumed to be normal with mean Z 
and coefficient of variation c,, where Z is the office average carried 
horizontal group load and where c, is the group-to-group coefficient 
of load variation for the office. The term c, may be inferred from office 
load balance data either by using analysis of variance or, more com- 
monly, by subtracting a standard value of the residual variance from 
the total measured variance of the office horizontal group load 
distribution. 


* This is the same value that was used in step-by-step dial-tone-delay calculations. 
(See Ref. 7.) These calculations are based on the Palm delay model® using an assumed 
mean call holding time of 150s and an assumed j factor of 5. In the notation of Ref. 
9, this value of 7 corresponds to a mean-time-to-defection of 30s. In Ref. 9, it is also 
stated that this value was found to be slightly conservative for most applications, 
based on review of panel office data reported in Ref. 9 and other (unpublished) 
step-by-step data. 
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To obtain the total offered blocked-call load for the office, we write 
the offered blocked-call load for an individual horizontal group x(z) 
as a function of the individual horizontal group carried load z and 
integrate the function over the office distribution of carried load. We 
now describe how 2(z) is calculated. 

Let N denote the number of subscriber lines per horizontal group. 
As discussed above, NV has been taken to be 35 in all computations. 
Let \, denote the combined originating and terminating rate when 
n subscriber lines are occupied and let uw, denote the subscriber line 
hang-up rate when n subscriber lines are occupied. Then 


An = A(N — n) Osns9 
=> (Nn) 10snsN-1 
and 
p= nA 0O<sn<10 
= 10/H + (n — 10)/A, lisnsvN, 


where H denotes the office average line-link holding time and H, de- 
notes the reciprocal of the defection rate for a call that is waiting for 
a line link to become available. As discussed above, H, is taken to be 
30 s based on results in Ref. 9. 

The parameter \ (combined originating and terminating rate per 
unoccupied subscriber line) is an unknown whose value will be ob- 
tained from the horizontal group-carried load z. The factor of 3 appear- 
ing in the definition of \, for 10 S$ n S N — 1 reflects the fact that, 
when all 10 line links are busy, an incoming call cannot cause a sub- 
scriber line to be occupied. The definition of the hang-up rate up for 
11 Sn SN reflects the assumption that, when all 10 line links are 
occupied, the holding times of the subscriber lines for which no line 
links are available should be shorter than full call holding times. 

For a horizontal group with carried load 2z, let 


7, = steady-state probability that n (out of NV) subscriber 
lines are occupied. 


Then, using standard methods for computing the steady-state dis- 
tribution of a birth-and-death process,” 


Nort *An—1 


) for = 1, 
Bist thn 


ro= cand my = c( 


where the constant c is determined by the requirement that 


N 
yay = bh 
n=0 
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Then the carried horizontal group load z is given by 


9 N 
= u nn + 10 2X Tn. (14) 
It is not hard to show that the carried horizontal group load z is a 
strictly increasing function of the subscriber line occupancy rate \ and 
that z < min [10, \NH]. This makes it easy to determine numerically 
the unique value of \ corresponding to a given carried load z. 

Once \ has been determined, the quantities An, un, and 7, are used 
to compute the offered blocked-call load contributed by the hori- 
zontal group. We take x(z) to be the expected number of occupied 
subscriber lines for which no line links are available. This yields the 
value 


x(z) = a (n — 10)rn. (15) 


IV. PEAKEDNESS OF THE BLOCKED-CALL STREAM 


In the model discussed in Section III, blocked calls arrive according 
to a Poisson process and defect after an exponential waiting time. 
Since blocked calls constitute an overflow stream and since it is well 
known!" that overflow traffic usually has a peakedness* greater than 1 
and hence is not Poisson, some discussion of the peakedness of the 
blocked-call stream is in order. 

The blocked-call stream is the superposition of overflow traffic from 
all of the (typically 400 to 600} horizontal groups in an office. In the 
case of an office with ideal load balance (identical horizontal group 
loads), all horizontal groups would have equal expected contributions 
to the blocked-call stream and standard limit theorems would suggest 
that the blocked-call stream should be approximately Poisson. 

Comparisons of calculated and observed dial-tone delays discussed 
in Section V (covering measured pTM occupancies up to about 0.84) 
indicate that, when the blocked-call stream is assumed to be Poisson, 
the calculated delays generally fall in the midrange of applicable data. 
However, there is a large variability in observed dial-tone delays 
measured in the same office under nearly identical levels of prm 
occupancy, second-failures-to-match, and percent all-ors-busy. The 
presence of this variability may be regarded as evidence that, in some 
busy hours, the blocked-call stream may be peaked in nature. Peaked- 
ness of the blocked-call stream is capable of accounting for substanti- 


* The peakedness of a stream of calls is, by definition, the variance-to-mean ratio 
of the number of busy servers when the stream is offered to an infinite group of 
servers with independent identically distributed exponential holding times. 
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ally higher calculated dial-tone delay at a given level of ptm occupancy, 
second-failures-to-match, and percent all-ors-busy than would be 
calculated under the assumption that the blocked-call stream is 
Poisson. 

One would expect the blocked-call stream to be peaked whenever 
most of the blocked calls are contributed by a small number of highly 
overloaded horizontal groups. Data on the horizontal group load dis- 
tributions during individual busy hours are not available for the test 
discussed in Section V and would be impractical to obtain on an on- 
going basis in any office. In the absence of data from which one could 
deduce directly the blocked-call peakedness, a treatment has been 
made using some results which R. I. Wilkinson obtained in the course 
of formulating his “equivalent random” method of characterizing 
overflow traffic." 

Wilkinson" assumes that traffic arrives and departs according to a 
birth-and-death process in which the arrival rate is increased when- 
ever the number of calls in the system exceeds a nominal number and is 
decreased whenever the number of calls in the system is less than this 
nominal number. The equilibrium distribution of the number of calls 
in the system resulting from these assumptions is shown to be negative 
binomial and, hence, this distribution is completely determined by its 
mean and variance (or by its mean and peakedness). Wilkinson then 
shows that negative binomial distributions rather closely approximate 
true overflow distributions under a number of different conditions. 

The effects of peakedness of the blocked-call stream on dial-tone 
delay were explored using the limiting analytical model of Section II 
by replacing the (truncated) Poisson distribution of blocked calls by 
a (truncated) negative binomial distribution of blocked calls. (Note 
that this does not require recomputing the conditional delay distribu- 
tions.) In this manner, dial-tone delay distributions were calculated 
and compared with the observed dial-tone-delay distributions dis- 
cussed in Section V, assuming peakedness values of 2 and 4. The re- 
sulting delay curves approximated the higher delays observed for 
given occupancy parameters. These calculations indicate that much 
of the observed variability in dial-tone delays under nearly identical 
load can be explained by peakedness of the blocked-call load. 

It is not difficult to modify the birth-and-death model discussed in 
Appendix B so as to accommodate peaked blocked-call input using 
Wilkinson’s method of representing this input.” 

Thus, when a negative binomial blocked-call distribution is used in 
the limiting model, the resulting distribution of the good-call dial- 
tone delay may be regarded as the limit of a sequence of good-call- 
delay distributions corresponding to models with blocked-call input 
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having the peaked form suggested by Wilkinson." The limit is taken 
as the bad-call-arrival rates approach zero and the mean bad-call- 
waiting times approach infinity with their products approaching posi- 
tive constants. Incorporation of the negative binomial distribution in 
the limiting model is accomplished simply by replacing the terms of 
the truncated Poisson distribution in eq. (138) with the corresponding 
terms of the negative binomial distribution, using any convenient 
truncation. 


V. COMPARISON OF CALCULATED AND OBSERVED DIAL-TONE 
DELAY DISTRIBUTIONS 


5.1 Summary 


Theoretical dial-tone-delay distributions, calculated using the ana- 
lytical model discussed in Section II, were compared with dial-tone- 
delay measurements made in a field test. The test was conducted in 
a No. 5 crossbar office with 60 LLFs, 4 pTMs, 68 dial-pulse (DP) ors, 
and 68 multifrequency (MF) ors. The data discussed in this section 
are from the time period February through April, 1974. 

The main conclusions of this study are that (z) the calculated delays 
generally fall in the midrange of applicable data, (77) there is a large 
variability in observed dial-tone delays measured under nearly 
identical levels of prm occupancy, second-failures-to-match, and 
percent of all-ors-busy, and (777) the field data show a clear increase 
in the ratio of waste ptm usage to total prm usage as frame load 
increases. The observed amount of increase in waste DTM usage agrees 
with theoretical predictions. 

In this section, the manner in which the dial-tone-delay measure- 
ments were taken is discussed and the method used to obtain the 
calculated delays is outlined. Next, some of the sources of variability 
in No. 5 crossbar dial-tone delays are identified and an explanation is 
given as to why a large variability in observed delays should be ex- 
pected. Two data plots are given indicating, respectively, the effects 
of frame load on waste pTM usage due to second-failures-to-match and 
the effects of frame load on incoming-first-failures-to-match. Finally, 
the conclusions of the study are discussed. 


5.2 Methods of obtaining the calculated and observed delay distributions 


Hourly dial-tone-delay measurements were made by placing ap- 
proximately 900 test calls per hour, using a standard 3-s dial-tone- 
delay testing machine which had been modified to record the pro- 
portion of test calls with delays exceeding X seconds, for X = 0.5, 
1.0, 1.5, 2.0, 2.5, and 3.0. 
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The observed actual DTM occupancy fz for a single hour is expressed 
in terms of the measured parameters by the formula: 











, _. (Total prm peg count) (0.015) + (Measured seconds of ptm usage) ; 
ed (4) (3600) 
(16) 


The term in eq. (16) involving the pTm peg count is to adjust for the 
seconds of ptm usage which pTM usage measuring devices do not 
record. 

The (adjusted) prm holding time for each hour was computed by 
A. R. Thorne from the prm peg count, the all-ors-busy peg count, 
the adjusted ptm usage [the latter of which comprises the numerator 
of eq. (16) ], and an additional “light-traffic adjustment” (used when- 
ever the observed actual ptm occupancy is below 0.80). Most of the 
adjusted ptm holding times are between 0.28 s and 0.31 s. A prm 
holding time of 0.30 s is assumed in the theoretical dial-tone-delay 
distributions discussed in this section. 

To compare predicted dial-tone delays with measured dial-tone 
delays, it was necessary to infer the observed increment in DTM occu- 
pancy due to prm second-failures-to-match (pTM2Frm). This incre- 
ment, denoted by A, is taken to be 


ree (Total prmM2rTm peg count) (HBC + 0.015) 


(4) (8600) » (17) 


where 


HBC = ptm holding time during a second-failure-to- 
match (seconds). 


A value of 0.40 s is used for HBC, based on data obtained during an 
earlier dial-tone-delay field test.° The observed yood-call DTM occupancy 
pu 1s defined to be 


Dial-tone-delay distributions were calculated, using the analytical 
model, for a range of values of actual prm occupancy p, and good-call 
DTM occupancy px, Where the parameters p, and p, are defined in 
Section II. The specific manner in which the distributions were calcu- 
lated is discussed below. The results of these calculations were tabu- 
lated into a set of dial-tone-delay distributions, indexed by the pairs 
(ps, pu). To compare the calculated and observed dial-tone delays, 
measured dial-tone-delay distributions from data-collection hours with 
similar values of 6, and é, were plotted on a graph along with one 
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theoretical dial-tone-delay distribution selected from the tabulation 
discussed above, such that ps  f» and pu & fu hold for the values 
pv and f, corresponding to the observed delays shown on the graph. 

These graphs are shown in Figs. 3, 4, and 5. Table I lists data per- 
taining to each of the observed delay curves on each of the graphs. 
At the top of each graph are listed the actual ptm occupancy p» and 
the good-call prm occupancy p, used for the theoretical dial-tone-delay 
curve (the solid line) on the graph. Also listed are the ranges of the 
fp, and fz values corresponding to the observed dial-tone-delay curves 
on the graph. The plotting symbols on the graphs indicate measured 
dial-tone delays. The dotted lines are smoothing curves fitted to the 
measured delays by the computer plotting routine used to draw the 
graphs. 

The predicted dial-tone-delay distributions were obtained in several 
steps. First, values of the good-call origination rate per LLF \; were 


LITY THAT DIAL—TONE DELAY tS AT LEAST X SECONDS 





SECONDS OF DIAL—TONE DELAY 


Fig. 3—Calculated and observed dial-tone delays. Actual occupancy in calcula- 
tions = 0.75 (data: 0.740 to 0.758). Good-call occupancy in calculations = 0.70 
(data: 0.684 to 0.697). 
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Fig. 4—Calculated and observed dial-tone delays. Actual occupancy in calcula- 
tions = 0.82 (data: 0.795 to 0.837). Good-call occupancy in calculations = 0.76 
(data: 0.744 to 0.782). 


computed so as to produce good-call prm occupancies of p, = 0.60, 
0.66, 0.68, 0.70, 0.72, 0.74, 0.76, 0.78, and 0.80. The values of 1 were 
obtained numerically from the formula py = \4NT/2p given in Sec- 
tion III. Note that p, the erlang C probability of all-ors-busy, is a 
function of \1, whereas NV and T are constants. For the office in which 
the test was conducted, N = 30 and, as discussed above, 7 = 0.30. 
In computing p, an average or holding time of 10.25 s was assumed 
based on data from the test. In all cases, the calculated values of p 
were greater than 0.99, so the all-ors-busy condition has a calculated 
probability of less than 0.01 under the conditions to which these 
distributions apply. (The observed fractions of all-ors-busy were 
also below 0.01 during most of the hours of the test. Hence, ors do 
not appear to have caused much of the dial-tone delay observed in 
the test.) 

Next, the conditional delay distributions corresponding to these 
parameters were computed using the analytical model. Using eq. (12), 
values of x (total erlangs of blocked calls) corresponding to a range 
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Fig. 5—Calculated and observed dial-tone delays. Actual occupancy in calcula- 
tions = 0.82 (data: 0.795 to 0.837). Good-call occupancy in calculations = 0.76 
(data: 0.744 to 0.782). (Note that a blocked-call peakedness of 4 is assumed in the 
calculations shown above. Figures 3 and 4 are based on an assumed blocked-call 
peakedness of 1.) 


of values of p, were computed for each fixed value of p,. For each such 
pair of p, and p,, the good-call dial-tone-delay distributions were then 
computed by (13) using the value of x, obtained as discussed above, 
and the conditional delay distributions corresponding to p,. The re- 
sult of these calculations is a set of good-call dial-tone-delay distribu- 
tions indexed by the pairs (pp, px). 

These good-call dial-tone-delay distributions include the indzrect 
effect of bad calls in that they reflect the increased ptm congestion 
produced by the bad calls. As discussed earlier, the direct effect of bad 
calls is expected to be small and therefore may be accounted for in a 
somewhat approximate manner. To reduce the number of variables 
that need to be considered, a single blocking probability is used in 
lieu of integrating the blocking probability formula over an assumed 
distribution of expected horizontal group loads. For the purpose of 
computing an average blocking probability ps, the average waiting- 
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Table | — Dial-tone-delay data used for comparison with calculated delays 


Fig. : 70 o No. pts Q * 
Date Time Actual ptm | Good-Call og cos/LLF | 0.5s*} 10s }1.5s/]2.0s]25s | 30s 

No. Occupancy | Occupancy Test Calls AORB 
3/4/74 1500-1600 74.0 69.7 900 0.38 N.A. 370 120 73 53 42 33 
3/25/74 | 1500-1600 74.5 68.4 900 0.30 1344 354 106 62 37 24 15 
3 4/9/74 1530-1630 75.5 69.4 900 0.00 1815 298 44 9 5 4 3 
2/6/74 1500-1600 75.8 69.7 900 0.01 1438 314 41 17 11 9 8 
2/22/74 | 1600-1700 74.7 69.4 900 0.06 1351 341 71 28 16 ll 7 
4 2/6/74 1600-1700 80.9 74.4 900 0.61 1480 461 188 122 90 76 59 
and 2/15/74 | 1600-1700 80.4 76.1 900 0.92 13890 430 175 121 83 62 49 
5 4/8/74 1530-1630 83.7 78.2 900 0.25 1384 426 126 53 25 14 10 


* For each hour listed, the number in this column is the observed number of dial-tone speed (pts) test calls whose dial-tone delay exceeds 
0.5 s. The columns labeled ‘1.0 s’’ through ‘8.0 s’’ are defined analogously. 


time-until-defection of a bad call (scwr) is assumed to be 10 s.* A 
consequence of this assumption is that the total bad-call origination 
rate corresponding to x total erlangs of bad calls on two pTMs is 7/10 
bad calls per second. The total origination rate of good calls on all 
30 LLFs is 30\;1. The average blocking probability PB corresponding to 
the occupancy pair (ps, pz) is then defined to be the ratio of the bad- 
call origination rate to the total origination rate of good and bad calls. 
Thus, 


x/10 


PB = 3/(10) + 30d)’ 


(19) 


where x and ), are the values corresponding to the occupancy pair 
(po, Pu): 

The conditional delay distribution of bad calls should be nearly 
exponential with mean delay uT/(10 — a), where ut is the average 
call holding time and a is the carried horizontal group load in erlangs.* 
For each occupancy pair (ps, pu), the theoretical dial-tone-delay dis- 
tribution for all calls is then given by 


P(D >t) = (1 — pB)P¢(D > t) + pB {exp |= eal , (20) 


where P¢(D > t) is the good-call dial-tone-delay distribution cor- 
responding to the occupancy pair (ps, pu). 

In the theoretical dial-tone-delay curves shown in Figs. 3, 4, and 5, 
a = 0 is used in eq. (20). The effect on P(D > 3 s) of setting a = 0, 
rather than using a more nearly correct value for each graph, is less 


*The value of scwT = 10 s used in calculating PB is consistent with the mean- 
time-to-defection H, = 30 s used in the horizontal group blocking calculations in 
Section III. The difference between the numerical values arises because these two 
parameters are defined differently. For the purpose of calculating ps, it is assumed 
that a bad call arrives, remains waiting for dial tone an exponential length of time 
with mean Bewr, and then defects. In the horizontal group blocking model discussed 
in Section III, a call which finds all 10 line links busy may either defect or may 
eventually obtain an idle line link. The call contributes to the bad-call load during 
the time that it remains waiting for one of the line links to become available. Thus, 
the quantity in the horizonial group blocking model corresponding most nearly to 
the parameter BcwT is the mean time until a call that finds all 10 line links busy 
either obtains a line link or defects. Based on delay calculations for the horizontal 
group blocking model discussed in Section III, and assuming a mean-time-to-defec- 
tion of 30 s for calls that do not obtain a line link and a mean line-link holding time 
of 150 s, the mean time until a call that finds all 10 line links busy either obtains 
a line link or else defects is calculated to be between about 10 and 12 s for horizontal 
group carried loads in the range of interest. Thus, it is reasonable to take BcwT = 10s. 

This expression for the conditional delay distribution of bad calls should over- 
estimate their delays somewhat because the expression does not account for finite- 
source effects. Since the fraction of calls experiencing these delays (i.e., the fraction 
of bad calls) is typically less than 0.01, the effect on 3.0-s dial-tone-delay probabilities 
of neglecting finite-source effects is typically less than 5 X 107-4. For this reason, 
these effects are neglected in eq. (20). 
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than 8 X 10~‘ in all cases. This error is in the opposite direction from 
the error, of comparable magnitude, resulting from not including 
finite-source effects in the bad-call delay distribution. 

Since the analytical model is a discrete time model with a time step 
of T seconds, P(D > kT) = P[D = (K+ 1)T] for k = 1, 2, ---. In 
Figs. 3, 4, and 5, 7 = 0.30. Because of the way that dial-tone-delay 
measurements are taken, a call which is recorded as having a delay 
greater than ¢ seconds may actually receive dial tone within a few 
milliseconds after time ¢t. Thus, in comparing the theoretical dial-tone 
delays with the observed dial-tone delays, the observed fraction of 
dial-tone delays greater than t seconds is taken as representing the 
observed fraction of dial-tone delays greater than or equal to ¢t seconds. 
The theoretical delay curve plotted is the curve P(D = kT) for 
k = 1, 2, ---, interpolated so as to produce a smooth curve. 


5.3 Sources of variation ir observed No. 5 crossbar dial-tone delays 


Dial-tone delay in No. 5 crossbar is influenced by a number of 
factors capable of producing a large variation in delays measured in 
different hours within the same office under very similar conditions 
of pT occupancy, percent all-ors-busy, and second-failures-to-match. 
As discussed in Section IV, much of the variability in dial-tone delay 
measured under very similar load conditions can be explained by 
differences in peakedness of the blocked call stream. Whenever most 
of the blocked-call load comes from a small number of extremely over- 
loaded horizontal groups, the blocked-call stream should have a peaked- 
ness greater than one. When a large number of moderately overloaded 
groups contribute to the blocked-call load, the blocked-call stream 
should be approximately Poisson (peakedness equal to one). Thus, 
differences in the individual busy-hour-load balance would be expected 
to produce different amounts of blocked-call peakedness, which in 
turn can account for appreciable differences in dial-tone delay mea- 
sured under nearly identical average load conditions. For example, 
Fig. 5 shows dial-tone delays calculated assuming a blocked-call peaked- 
ness of 4 for the same conditions shown in Fig. 4, which is based on a 
blocked-call peakedness of 1.* The calculated delay curve in Fig. 5 fits 
the top two observed delay distributions rather closely. 

Some additional identifiable sources of dial-tone-delay variation 
under similar load conditions are: (7) within-hour trends in traffic, 
(22) nonstandard (and possibly erratic) gating caused by improper 
functioning of the master traffic controller circuitry, (i) ptm prefer- 


“The method by which blocked-call peakedness is treated in the model is dis- 
cussed in Section IV. 
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ence for calls from a small subset of lines on each horizontal group, 
(wv) variation in ptm first-failures-to-match, and (v) competition be- 
tween DTMs and completing markers for line-link connectors. 

The first of these sources should produce effects similar to those of 
‘blocked-call peakedness. The second source may cause nonuniform 
congestion. The third source is predicted to result in a slight outward 
shift in the delay curve. The fourth and fifth sources should be reflected 
in increased measured ptm holding time and in increased pTM occu- 
pancy. Although approximate allowances can be made for the average 
congestion increase produced by some of these phenomena, no quanti- 
tative estimate is available for their total contribution to hourly 
variation in dial-tone delay. 

In addition to identifiable sources of dial-tone-delay variation, simu- 
lation studies indicate that there can be an appreciable residual varia- 
tion in simulated hourly No. 5 crossbar dial-tone delays obtained in 
different runs with identical input parameters (and, hence, with identi- 
cal expected load conditions).* 

Figure 6 shows four dial-tone-delay distributions obtained using 
the gating simulation model. In this model, the blocked-call stream 
is Poisson. These distributions were produced by simulating four 
individual hours, using identical input parameters. The delay distribu- 
tions shown are for the calls that did not encounter horizontal group 
blocking and are based on the total number of such calls processed 
during the hour. The set of four 3-s dial-tone delays has a coefficient 
of variation of 0.28 and a mean of 0.069. (The coefficient of variation 
is the ratio of the standard deviation to the mean.) Plotted on the 
graph along with the delay curves are error bars indicating the 2-sigma 
limits of 0.034 and 0.107 associated with the above mean and coeffi- 
cient of variation. 

Actual dial-tone-delay measurements are based on test calls. During 
a given busy hour in a typical No. 5 crossbar office, approximately 
900 test calls are made on a fixed set of 60 (out of 600) horizontal 
groups. The use of test calls introduces sampling error, which is not 
represented in the distributions shown in Fig. 6 and which would 
have an associated coefficient of variation of about 0.12 for the parame- 
ters applicable to Fig. 6. 


5.4 Data on frame load effects 


Figure 7 is a data plot of line-link frame load versus waste pT usage 
due to second-failures-to-match based on data from the test. The 
quantity “‘prm usage fraction due to 2FrMs’”’ shown on the ordinate is 


*-This result is one of the main conclusions of an earlier No. 5 crossbar dial-tone- 
delay simulation study conducted by S. Halfin? 
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PROBABILITY THAT DIAL—TONE DELAY EXCEEDS X SECONDS 





0 1 2 3 4 5 6 
SECONDS OF DIAL—TONE DELAY 


Fig. 6—Simulated No. 5 crossbar dial-tone delays. Distributions are based on simu- 
lation of four individual hours with identical inputs. 


defined by 
Total ptm usage due to 2FTMs 
Total prm usage 





DTM usage fraction due to 2FTMs = 


Pile 


Figure 8 is a data plot of line-link frame load versus incoming-first- 
failure-to-match (1FFM) based on data from the test. In each figure, 
each of the data points represents data from 1 h. The average actual 
DTM occupancy (averaged over all the data-collection hours with 
measured line-link frame loads of 1100 ccs/LLF or more) is 0.54; the 
DTM occupancy range is from less than 0.40 to 0.84. 

Existing theory indicates that, for a given office configuration (in- 
cluding a given junctor pattern), 1rrM is directly and primarily de- 
pendent on frame load. This conclusion is borne out by Fig. 8, which 
demonstrates a well-defined trend (with a moderate amount of data 
scatter) of increasing IFFM with increasing frame load. 

The extent to which frame load affects prm usage and dial-tone 
delay in any given hour depends, in an indirect way, on several differ- 


DIAL-TONE DELAY 1783 


DTM USAGE FRACTION DUE TO 2FTMS 





0 
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CCS/LLF 


Fig. 7—Effect of frame load on increased ptm usage due to 2FTMs. The assumed 
pT holding time during a 2rTm = 0.40 s. 


ent variables, including ptm occupancy and the distribution of carried 
load among horizontal groups. Because of this dependence, the nature 
of which has been discussed more fully in previous sections, the data 
scatter in the plot of waste ptm usage fraction versus frame load (Fig. 
7) is much larger than in the plot of 1rrm versus frame load (Fig. 8). 

The curve labeled ‘THEORETICAL’ in Fig. 7 was calculated using the 
limiting analytical model [eq. (12) ] and the horizontal group blocking 
model discussed in Section III. In these calculations, the blocked-call 
load was assumed to be Poisson (peakedness equal to 1). 

To obtain the theoretical curve, it was first necessary to determine 
what calling rates should be assumed in the calculations. The calling 
rates were inferred from the data upon which Fig. 7 is based by first 
using linear regression (least squares) to express the observed good-call 
DTM occupancy, 6~, aS an empirical function of the observed frame load, 
ccs/LLF. For each value of the frame load, the calling rate was taken 
to be the particular calling rate corresponding to the least squares 
value of #1, assuming an average pTm holding time of 0.30 s and a 
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INCOMING FIRST—FAILURE—TO—MATCH 
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Fig. 8—Iiffect of frame load on incoming first-failure-to-match. 


zero probability of all-ors-busy. The justification for the assumptions 
regarding the prm holding time and the probability of all-ors-busy 
is given earlier in this section. 

Next, the horizontal group blocking model discussed in Section III 
was used to compute the expected blocked-call offered load correspond- 
ing to each value of frame load and (empirically associated) calling 
rate. The office distribution of carried horizontal group load was 
assumed to be normal with group-to-group coefficient of variation 
inferred from the office horizontal group load distribution using the 
method discussed earlier. The calculated values of offered blocked- 
call load were then used in eq. (12) to compute the theoretically pre- 
dicted fraction of waste DTM use, (ps — px)/ps, corresponding to these 
frame loads and calling rates. 

Since both the theoretical and the observed values of waste DTM 
use depend not only on frame load but also on the calling rate, neither 
the data plotted in Fig. 7 nor the theoretical curve shown on the figure 
should be regarded as being applicable to calling rates or load balance 
conditions other than those upon which this figure is based. 
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The 


purpose of Fig. 7 is to illustrate—for the conditions of calling 


rate, frame load, and load balance represented in this study—an 
empirical relationship between frame load and the ratio of waste pTmM 
usage to total ptm usage with increasing frame load and to show that 
this empirical relationship is consistent with predictions of the theo- 


retical 


models. 


5.5 Conclusions 


The 


conclusions which the author has drawn from this comparison 


of theoretical and observed dial-tone-delay distributions are as follows: 


(2) 


(it) 


(222) 


(wv) 
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Where two or more observed delay distributions appear on the 
same graph, the theoretical delay distributions usually fall 
approximately midway between the maximum and minimum 
observed delay distributions. This indicates that, in the 
(actual) pTM occupancy range spanned by these data (pTM 
occupancies up to about 0.84), the analytical model shows 
good agreement with the data. 

Much of the large variability in observed 3-s dial-tone delays 
measured under nearly equal conditions of pTM occupancy, 
second-failures-to-match, and percent al]l-ors-busy can be ex- 
plained by assuming the blocked-call stream to have different 
peakedness values (ranging from 1 to about 4) in different 
hours. High blocked-call peakedness, illustrated by Fig. 5 in 
which the peakedness is taken to be four, would be expected 
whenever most of the blocked-call load comes from a small 
number of extremely overloaded horizontal groups. Low 
blocked-call peakedness, illustrated by Fig. 4 in which the 
peakedness is taken to be 1, would be expected whenever a 
large number of moderately overloaded horizontal groups con- 
tribute more or less equally to the blocked-call load. 
Simulation results indicate that there should be a large residual 
variability in 3-s dial-tone delays measured in different 
(simulated) hours under identical expected load conditions. 
This variability, illustrated in Fig. 6, is in addition to the 
variability due to blocked-call peakedness discussed above. 
A third nonnegligible source of variability in observed dial- 
tone delays is the sampling variability associated with the use 
of test calls to measure delays. 

For the conditions of calling rate, frame load, and load balance 
represented in this study, the observed increase in the ratio of 
waste DTM usage to total ptm usage as frame load increases is 
consistent with theoretical predictions. This is illustrated in 
Fig. 7 and discussed in Section 5.4. 
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APPENDIX A 
Computation of the Two-Marker Stationary Distribution for k = 1 


The purpose of this appendix is to prove that the. limit r defined 
by (10) exists and that the quantities v; = a; + 78;, satisfy (21), 
where a; and 8; are defined by (9). 


0; = 0, all 7 =0 
Vo = 1; 
ite 
v= LvP.(i, 7), allj 20 
i=0 
x vy <, (21) 
j=0 


We first show that there exists exactly one positive number r such 
that v; = a; + r8; satisfies (21). Then we show that the limit in (10) 
exists and is equal to this number. 

When \ < 2p, we know that the stationary distribution P; exists 
and P;(j) > O for each j. Let 7; = P,(j)/Px(0). Then for r = P;(1)/ 
P,,(0), we have 

To = dy 


71 = 7, 


J+2 
wr; = > Pili, j), all f= 0. 
i=0 : 


Now let v; = a; + "Bi, with r as defined above. Then, it follows from 
(9) that 


Vo = 1, 
y= 7, 
j+2 
= >) oP, (2,7), allj 2 0. (22) 
ix0 


Since vo = mo and vy = 7m; and since the vectors v and = both satisfy 


DIAL-TONE DELAY 1787 


the same recurrence relation, which has the property that the jth 
term for each j = 2 is uniquely determined by the Oth and Ist terms, 
we see that v; = 2; > 0 for all 7 = 0. Since the Markov chain is posi- 
tive recurrent for \/2p < 1, it follows from (22) and Karlin™ that 
> j=0 0; <«. Hence, v satisfies (21). Thus, there exists at least one 
value of r for which v; = a; + r§; satisfies (21). If r’ is any number 
such that v; = a; + 7’8; satisfies (21), then by uniqueness of the 
stationary distribution, 7’ = v;/vp = P.(1)/P.(0), so r’ = r. Thus, 
exactly one such r exists. 

Since a; + r8; > 0 for each 7, it follows that (10a) holds. Hence, 
the increasing sequence m, is bounded above by r and the decreasing 
sequence JZ, is bounded below by r. So both m = lim, m, and 
M = lim, M, exist and satisfy 


msrsnM, (23) 
Mn SM, 
Ms M,. 4) 


If either of the inequalities in (23) were strict, then m < M and, in 
view of (24), each x in the interval m < x < M would satisfy 


Mm S2sM, 


for each n, from which it follows that v; = a; + x6; satisfies (21). 
This contradicts the uniqueness of r and proves that m = M = r. 

It is easy to extend the above result so as to give a method for 
computing the stationary distribution of any irreducible positive re- 
current Markov chain on the nonnegative integers such that a fixed 
positive integer no exists for which P(j + no, 7) > 0 and P(j +n, J) 
=0 for all 7 and all n > no. 


APPENDIX B 
Convergence in Distribution to the Limiting Model 


In the limiting model, the number of bad (i.e., blocked) calls in the 
system is a time-independent, truncated-Poisson-distributed random 
variable. We will show that the steady-state distributions of good- and 
bad-call queue lengths and the steady-state good-call dial-tone-delay 
distribution in the limiting model are the limits of the corresponding 
distributions for a sequence of models in which the bad call queue 
lengths form birth-and-death processes. 

In the nth model, the bad-call arrival rate, denoted by A‘, ap- 
proaches zero and the bad-call mean waiting time until defection, de- 
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noted by H‘, approaches infinity in such a way that their product 
(total expected erlangs of bad calls) approaches a finite positive con- 
stant x. 

The physical motivation for considering the limiting model is that 
the actual situation is one in which blocked calls appear infrequently 
(relative to total call arrival rates) and tend to remain in the system 
for a long time (relative to ptm holding times). The mathematical 
motivation is that the limiting model is much easier to analyze than 
a model in which the bad-call arrival and departure processes are 
represented explicitly. 

Preliminary computations using a single marker version of the 
analytical model, which represented the bad-call arrival and departure 
processes explicitly, showed that different bad-call arrival rates and 
waiting times had little effect on dial-tone-delay distributions as long 
as the product of the bad-call arrival rate and bad-call waiting time 
remained constant. In addition, these distributions were all quite 
close to those obtained from the limiting form of the model. These 
results make sense intuitively because bad calls simply cycle through 
the system, absorbing some ptm uses while present, and eventually 
defect; thus, what should matter is mainly the distribution of the 
number of bad calls in the system at any time. In the version of the 
analytical model discussed in this appendix, the blocked call queue 
length distribution is shown to be a truncated Poisson with mean 
equal to the total erlangs of bad calls. 

Let K be the maximum possible number of bad calls that can be 
in the system at any time. In all computations using realistic No. 5 
crossbar busy-hour input parameters, the expected erlangs of bad 
calls have been low enough that the Poisson probability of more than 
five blocked calls being present in a two-marker system has been less 
than 10~%. Hence, K may be taken to be 5. (Note that in an actual 
No. 5 crossbar office, the maximum number of bad calls that can be 
in the system at any time is trivially bounded above by the number of 
subscriber line terminations in the office.) 

Let Sn(p) denote the saturation load of the nth model. Since no 
more than K blocked calls can be in the queue at any time, in the 
nth model, we have 2p = S2(p) = Sn(p) = S2(p) = 2p for all n; 
hence, S,(p) = 2p. Thus, the saturation load of the nth model is the 
same as the saturation load of the limiting model. For the nth model, 
let 


I 


X*? = number of good calls in the queue at time mT, 
Y() (mT) = number of bad calls in the queue at time mT, 
pilk, lt) = P{YME(m + OT] = UY (mT) = k}. 
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In the nth model, Y™(é), the bad-call population size at time /, is a 
finite-state space birth-and-death process with birth-and-death rates 


he = Am O<msK-—-1 
= otherwise 

OsmskKk 

= otherwise. 


Since the queuing model is in discrete time, only the values of Y“ (¢) 
at t = mT are of interest. The transition matrix p(k, 1) is given by 


— TA 
p=e ", 


where A, is the infinitesimal generator matrix of the birth-and-death 
process Y‘™, (See Ref. 15.) Since lim,..,, \™ = Oandlim,.,, H™= 0, 
it follows that lim,.,. An = 0. Hence, 


lim pa(k, tl) = 8x1 (25) 
for each k, 1 = 0,1, ---, K, where 6,. = O fork ¥ land 6, = 1. 
We also have 


P{XM. = 97, YRL(m + YT] = XP = 1, VRr = k} 
= Prlk, lI) Px(t, D: 
where P;(2, j) is given by (1) and (2). Since the nth model has satura- 


tion load 2p, we know that a stationary queue length distribution 
an(j, l) exists for \/2p < 1 and satisfies the equation 


ra(D = XX walé, Dp alb, )Pulé, 5). (26) 


Let 7,(-, 1) be the marginal equilibrium distribution of the number of 
bad calls in the system. Then 


ao) =X wal, D (27) 
and by (26) 
a(-D = 2 wali, Hpalk, DE Pali, 1) 
= ¥ wal-, Walk, 0, (28) 


Hence, 7n(-, k) is the stationary distribution of a Markov chain whose 


1790 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1975 


transition matrix is p,(k, l). Using standard theorems on birth-and- 
death processes (Ref. 10) it follows that, for all values of / = 0,1,---, K, 


ri(-,k) = lim PLY™(t) = k/Y™(0) = 7] 








~ CraKtn : 
aes © O<kSK 
= 0 otherwise, (29) 
where 
Ln = MYH(™) 
and 
K gt 
CronK = kl 
We will now show that 
: ; Crpx* 


for allk = 0,1,---, K andall 7 = 0,1, ---, where P; is the stationary 
distribution of good-call queue length in the model with k bad calls 
permanently present in the system. The right-hand side of (30) is the 
queue-length distribution for good and bad calls in the limiting model. 

Let X‘) denote the number of good calls in the queue for the nth 
model in equilibrium. Since the number of bad calls present in the 
nth model is always less than or equal to K for all n, it is clear on 
intuitive grounds and can be proved rigorously using stochastic order- 
ing that 


PIX™ 2 j] SY Pr) (31) 


for each j. 
To prove (80), it suffices to show that if 7, is any subsequence of 
awn for which lim, tn/(j, k) = (J, &) exists for each 7 and k, then 
it 


(i) (5, b) = ¥ w(G, E)Pr(i, 3) 


1=0 


and 





2 —1,.k 
(i) w(j,k) = SE. 
fo kl 


To see that (z) and (72) are sufficient, note that, by uniqueness of 
the stationary distribution, P;, (¢) and (iz) imply 


Crna 


r(j,k) = GEE P,(j). (32) 
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Since 0 S 7,(j, k) S 1 for all j and k, every subsequence 7, contains 
a further subsequence 7," for which the limit lim, 7,//(j, k) exists 
for all j and k. In view of the above discussion, all of these subse- 
quences have the same limit, namely the right-hand side of (32). 
Equation (80) follows. 

We now prove (2) and (zz). Condition (2) follows immediately from 
(25) and (26). To see that (7) holds, we proceed in two steps. First, 
note that by (27), (29), and Fatou’s lemma, 


fe) 


» (4) S < CeKe 





for each k. If any of the above inequalities were strict, then 


Me- 
Ms 


ll 
° 


(I; k) <1. 


k 


Hence, to prove (72) it suffices to show 


i 4 


» (i, k) 2 


ie., we must show that, in the limit, no probability mass escapes to 
infinity. Let « > 0 and choose jo such that 


¥ Px(j) <« 
Then, using eq. (31), 


= lim 2 53) Tn (2, k) 


nit og 


= lim P[IX@") < jo] 


jo 
= X& Px) 
Zl-—e. 


Thus, (77) holds and eq. (30) follows. 

It can be readily shown, using the above results together with 
standard theorems, that the delay probabilities in the nth model 
converge to those in the limiting model. (See Refs. 16 and 17.) 
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